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Abstract 

We show that many observations of W44, a supernova remnant in the galactic plane at a 
distance of about 2500 pc, are remarkably consistent with the simplest realistic model. The 
^ - model remnant is evolving in a smooth ambient medium of fairly high density, about 6 cm~^ 

on average, with a substantial density gradient. At the observed time it has an age of about 
20,000 years, consistent with the age of the associated pulsar, and a radius of 11 to 13 pc. 
Over most of the outer surface, radiative cooling has become important in the post shock gas; 
' on the denser end there has been sufficient compression of the cooled gas to develop a very 

QQ I thin dense half shell of about 450 Mq , supported against further compression by nonthermal 

. pressure. The half shell has an expansion velocity of about 150 km s~^, and is bounded on 

the outer surface by a radiative shock with that speed. 
— ,: , The deep interior of the remnant has a substantial and fairly uniform pressure, as expected 

I ■ from even highly idealized adiabatic models; our model, however, is never adiabatic. Thermal 

O ' conduction, while the remnant is young and hot, reduces the need for expansion cooling, and 

^ , prevents formation of the intensely vacuous cavity characteristic of adiabatic evolution. It 

^2 ' radically alters the interior structure from what one might expect from familiarity with the 

^ ' Sedov solution. At the time of observation, the temperature in the center is about 6x10^ 

K, the density about 1 cm~^. The temperature decreases gradually away from the center, 
1 J - while the density rises. Farther out where cooling is becoming important, the pressure drops 

r> I precipitously and the temperature in the denser gas there is quite low. Our model is similar 

^ . to but more comprehensive than the recent one by Harms et al. (1997). Because their model 

lacked thermal conduction, ours is more successful in providing the thermal x-rays from the 
hot interior, including a better match to the spectrum, but neither provides the sharpness of 
the central peaking without further complications. 

By using a 2d hydrocode to follow the evolution in a density gradient, we are able to verify 
that the spatial and velocity structure of the HI shell are a good match to the observations, 
without the complications suggested by Koo and Heiles (1995), and to demonstrate that the 
remnant's asymmetry does not substantially affect the distribution of x-ray emitting material. 
A Id hydrocode model is then used to explore the effects of nonequilibrium ionization on the 
x-ray spectrum and intensity. We calculate the radio continuum emission expected from the 
compression of the ambient magnetic field and cosmic rays into the dense shell (the van der 
Laan mechanism, 1962a) and find it to be roughly consistent with observation, though the 
required density of ambient cosmic ray electrons is about 4 times greater than that estimated 
for the solar neighborhood. We estimate the optical emission that should be present from 
fluorescence of UV, emitted by the forming shell and the radiative shock and absorbed in 
the cold shell and the ambient medium, and the associated 63 fim [OI] emission. Both are 
in agreement with the intensity and spatial structures found in recent observations. Neither 
requires interaction with a dense molecular cloud for its generation. We calculate the gamma 
rays that should be emitted by cosmic ray electrons and ions in the shell, interacting with the 
cold material, and find each capable of generating about 25% of the flux reported by EGRET 
for the vicinity. 

We provide several analytic tools for the assembly of models of this type. We review 
the early evolution and shell formation analyses and their generalizations to evolution in a 
density gradient. We also calculate the density and temperature that should be present in the 



hot interior of a remnant with thermal conduction. We supply the van der Laan mechanism 
in a particularly useful form for the calculation of radio continuum from radiative remnants. 
Finally, we demonstrate a simple technique for estimating the optical emission expected. These 
tools are employed to choose parameters of models which we then explore with our Id and 2d 
hydrocodes, providing, respectively, the detailed x-ray spectra and dynamical characteristics. 
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— masers — radio continuum — interstellar matter 
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1 Introduction 

W44 (3C392) is a supernova remnant (SNR), located near the Galactic plane {I = 34.7° , b — 
—0.4°), which has a centrally peaked distribution of thermal x-ray surface brightness, coupled 
with a shell-like radio continuum morphology. This combination makes it a member of the SNR 
class which also includes W28, 3C400.2, Kes 27, MSH 11-61A, 3C391, and CTBl (Rho 1995). 

In this paper we demonstrate that many observed features of this remnant are consistent with 
the simplest realistic model, a remnant evolving in a fairly smooth ambient medium of moderate 
density with a substantial density gradient, and with thermal conduction active in its interior. It 
has a radiative shock and dense shell on one side which is emitting copious amounts of optical and IR 
radiation, and is sweeping up and compressing ambient cosmic rays and magnetic field to produce 
radio synchrotron and gamma rays. This interpretation differs substantially from earlier attempts 
to model the remnant's x-ray emission or 21 cm characteristics as arising from a population of 
evaporating clouds in the interior or the disruption and acceleration of a wind blown shell, but is 
similar in spirit to the recent investigation by Harrus et al. (1997). Our model differs from the 
latter in that it includes thermal conduction within the hot interior of the remnant, assesses the 
effects of the density gradient, and is considerably more detailed dynamically and more ambitious 
in comparison with observations. Our interpretation also differs from past efforts to understand 
the radio continuum, the infrared emission, and associated OH masers: we find no evidence for 
interaction of the remnant with dense molecular clouds in the vicinity. 

Section ^ opens with a discussion of the complexity of the region in which W44 is located, 
followed by a revision in the estimated distance (using a recent rotation curve and the standard 
galactocentric distance of the Sun). It then surveys the observations of W44, including radio 
continuum, 21 cm. X-rays, characteristics of the pulsar and its weak nonthermal nebula, optical. 



IRAS and ISO observations, magnetic field measurements via polarization of OH masers, and the 
gamma rays seen with EGRET. 

Section ^ supplies the analytic tools for estimating remnant parameters consistent with the 
observations. It reviews the Sedov and cooling phases for uniform density, provides a means for 
estimating the effects of thermal conduction on the characteristics of the remnant interior, surveys 
the anticipated shell characteristics, introduces the radio continuum expected from the van der Laan 
mechanism, and provides simple tools for estimating the Ha, 63 /zm, and gamma ray intensities. It 
then shows how, using the approximations of Maciejewski and Cox (1998), the remnant parameters 
can be estimated in the presence of a density gradient, closing with a summary of what can be 
learned from analytical modeling alone. 

Section ^ presents the physical assumptions of the hydrocode modeling and describes their 
implementation and justification. This includes the use of a nonthermal pressure term, the thermal 
radiation emissivity and spectrum, and the assumptions regarding dust, thermal conduction, and 
synchrotron emission. 

Section ^ presents the results of 2d modeling, exploring the density, temperature and pressure 
structures of the hot interior and the mass, velocity, and nonthermal emissivity structures of the 
dense shell. Extensive comparison is made with both the 21 cm and radio continuum observations. 
The non-equilibrium ionization could not be followed, but we verified that the structure in the hot 
interior is nearly spherically symmetric and should be reasonably approximated by the Id model, 
which is able to track the ionization of the plasma. 

Section ^ presents the results of Id modeling, aimed principally at obtaining the x-ray surface 
brightness distribution and its spectrum, and comparison with those observations. 

Finally, Section ^ provides the discussion and conclusions. 

2 Survey of Observations 

The essential observed properties of W44 itself are summarized in Table |l] at the end of this 
section. Detailed description beyond the survey below can be found in recent papers by Giacani 
et al. (1997), de Jager & Mastichiadis (1997) and Harms et al. (1997). 

2.1 Environment 

W44 lies in a complex region of the sky, in close proximity to several molecular clouds (Sato, 1986), 
in a giant molecular cloud (GMC) complex (Dame 1983, Dame et al. 1986), and at the base of the 
Aquila Supershell (Maciejewski et al. 1996). The ambient density in which the remnant evolves 
is fairly high (averaging about 6 cm^"^, from our results below), presumably because of its GMC 
environment. Filamentation of the surface emission shows that the medium is moderately irregular 
on scales somewhat smaller than the remnant. It has been suggested that the remnant may be 
impacting directly on a very dense molecular cloud; the interpretation of the clearly associated 
OH masers by Claussen, et al. (1997), for example, is based on that assumption. Our model does 
not include this possible interaction and appears not to require it to understand any of the other 
observations. Considering the superbubble and its funnel-like connection to the GMC complex in 
the galactic plane, it is tempting to suppose that W44 might be the latest of a series of SNe in 
the general region, the earlier ones having contributed energy to the superbubble. Thus far, there 
is no evidence for any transport of energy from W44 except by cooling radiation and, possibly, 
escaping cosmic rays; its high density environment has quite effectively localized the mechanical 
and thermal components of its energy. Any contribution it may make to the superbubble lies in 
the future. 

2.2 Distance Determination 

In order to estimate the distance to the supernova remnant, Caswell et al. 1975 measured the 
21 cm absorption profile toward W44. They found components at velocities of -1-12, -1-22, -1-30, 
and -f42 km sec^^. (In confirmation, Goss, Caswell, and Robinson (1971) observed 1667 MHz 
OH absorption components at -1-12.5, -1-22, -1-30, and -1-42 km sec~^.) Because there are no higher 



velocity components, the +42 km sec~^ component corresponds to the near kinematic distance. 
Using the Schmidt (1965) galactic rotation curve and 10 kpc as the distance between the Earth and 
the galactic center, Caswell et al. estimated the distance to W44 as 3 kpc. We have recalculated 
the distance using the Clemens (1985) galactic rotation curve and 8.5 kpc as the distance between 
the Earth and the galactic center. The new estimate is 2.5 to 2.6 kpc. 



2.3 Radio Continuum Observations 

W44 appears as an elongated shell-type remnant in radio continuum maps. The size, measured 
from the Giacani et al. (1997) 1442.5 MHz image is 25 x 35 . In the less detailed map from the 11 
cm survey of Reich et al. (1990), the remnant appears noticeably brighter over the eastern half. 

The radio continuum image presented by Jones, Smith, & Angellini (1993) shows a very highly 
filamentary structure indicative of a severely edge brightened distribution of emission. The fila- 
ments are concentrated on the eastern side and are contained within a sharp boundary, while the 
western side contains only one very bright feature and has a muted boundary. 

Kundu and Velusamy (1972) show that W44 has a very regular linear polarization pattern. 
In the brighter half of the source, the polarization reaches 20%, and the magnetic field vector is 
almost parallel to the long axis of the remnant. The bright filaments tend to lie along this axis as 
well. 

The flux densities have been measured over a wide range of radio frequencies (e.g. by Kundu and 
Velusamy, 1972; Dickel and DeNoyer, 1975; Clark and Caswell, 1976; Clark, Green, and Caswell, 
1975; Kassim, 1992). The spectral index is about -0.33 (Kovalenko et al. 1994) to -0.4 (Giacani, 
et al. 1997), the latter reporting a flux density of about 400 Janskys at 330 MHz. Table |l| shows 
500 Janskys at 100 MHz as representative. 



2.4 Hi 21 cm Emission from W44 

Koo & Heiles (1995) have made a 21 cm velocity resolved map of W44 and found emission up to 
i^LSR = 210 km s^^. By subtracting off the systemic velocity of the clouds in this region (~43 
km/s, Sato 1986), they were able to estimate the shell's expansion velocity as 150 ± 15 km s~^. 
The observations at low velocity near the limb were badly confused with ambient material, leaving 
an accessible velocity range only above 130 km s^^ LSR. The atomic hydrogen mass in the shell 
integrated over this range was about 73 Mq . Thus, if the neutral material is confined to a thin shell, 
the average HI column density of the receding side of the remnant is about 3 x 10^^ atoms cm~^. 
Koo and Heiles searched for the near side of the SNR, but found no emission with wlsr < — 70 km 
s-i. 

By fitting the profile of an expanding shell to their position-velocity diagrams, Koo and Heiles 
(1995) estimated the full angular size of the HI distribution from their observations over a restricted 
velocity range. They concluded that the HI distribution is somewhat smaller than and interior to 



the radio continuum shell. This conclusion will be disputed in section 5.2.2 below where the spatial 



distribution and position- velocity diagrams of the HI data are compared directly with our model. 



2.5 X-ray Observations 

The ROSAT PSPC (Rho et al, 1994) and Einstein IPC (Smith et al, 1985) x-ray maps show 
centrally concentrated thermal emission whose dim outer edge extends as far as the radio continuum 
shell in the north, south, and west (Rho et al. 1994). (This appearance, however, may be due in 
part to scattering of the x-rays by the substantial column density of dust along the line of sight.) 
The X-ray emitting region may similarly extend to the eastern edge where it is absorbed by a 
foreground molecular cloud which covers irregularly from a(1950) = 18''53'"40* to 18''54™55^ and 
^(1950) = 01°10' to 01°20' (Wootten, 1977). Rho et al (1994) separately analyzed 11 spatial 
zones in the ROSAT PSPC image. From fitting the pulse height distribution with equilibrium 
ionization thermal models, they found values for the temperature of the X-ray emitting gas from 
3.9 to 7.6x10^ K. The fitted values for the column density of intervening absorbers ranged from 

1.6 to 2.1 xlO^^ cm~^. From the spatial distribution of these results they conclude that there is no 
significant radial trend in temperature. Working from the fitted results is somewhat misleading. 



however, as Figure 7 of Rho et al. (1994) clearly shows that the spectral hardness is significantly 
greater in the center where the x-rays are brightest. 

Smith ei aL (1985), Jones ei aL (1993), and Rho et a/. (1994) calculated the x-ray luminosity for 
the remnant interior from the Einstein, EXOSAT, and ROSAT data, respectively. Their estimates 
should be reduced by about 25 to 30% at our revised distance. The implied emission measure 
within the remnant requires an rms electron density of about 0.4 to 0.75 cm~^. 

The details of the observed x-ray surface brightness distribution and spectrum are presented in 
section |6| where they are compared with the results of our Id model. 



2.6 The Pulsar and its Synchrotron Nebula 

An apparently associated 267 millisecond pulsar (PSR B1853+01) lies projected on the southern 
part of the remnant (Wolszczan, Cordes, and Dewey 1991). Its age (^^20,000 years ±20%) provides 
an estimate for that of W44, and is typical of radiative remnants. The pulsar powers a small 
synchrotron nebula, observed in 4.0 to 9.5 keV ASCA data (Harms, Hughes, and Helfand, 1996), 
and in 4860 and 8440 MHz VLA data (Frail, Giacani, Goss, and Dubner, 1996). The x-ray flux 
from synchrotron is much smaller than the thermal component. 



2.7 Optical 

Giacani et al. (1997) have taken Ha and [SII] images of W44. All observed optical emission is 
within the boundaries of the radio shell, but at the sensitivity of the observation is not found 
everywhere within it. There appears to be noticeably higher foreground extinction in the south. 
The brightest Ha features have a surface brightness ~ 1 x 10~^^ ergs cm~^ s~^ arcsec"^ while 
the faintest features visible in the image have about a tenth that. The brightness of the [SII] lines 
is not reported but the exposure time was the same and the image showed comparable structure 
leading us to believe that the intensity is similar. Much of the optical emission is filamentary and 
in excellent correlation with radio continuum filaments. In some places there is fairly bright Ha 
that is smoother and not well correlated with the radio. On the eastern half of the remnant, where 
the radio is generally brighter, the observable optical emission is much lower. The Ha features 
correlate with [SII] emission, which often traces radiative shocks. Conversely, the spatial correlation 
between the optical and x-ray emission is poor, even on very large scales, differing entirely from 
what would be expected if the Ha and [SII] emission arose from evaporating clouds which were 
loading the interior with mass and enhancing the density of hot gas. 



2.8 Infrared, OH Masers, and Gamma Rays 

Giacani et al. (1997) show IRAS observations in 60 and 100 /im in a large field around W44, which 
they say indicate emission compatible with an origin in shock heated dust, while ISO observations 
presented by Reach and Rho (1996) indicate that, in all fields they examined, the continuum is 
characteristic of dust heated by the diffuse interstellar radiation field and is dominated by pre-shock 
and unrelated clouds. The IRAS map, however, shows weaker emission interior to the northern 
and mid-southern sections of the remnant not explored by ISO, that may well arise from hot dust 
within the remnant. (The ISO sampling crossed the southern edge of the northern enhancement, 
and Reach and Rho report that the continuum is brighter there and peaked at the edge where 
there is a radio ridge, but make no further comment about the heating mechanism.) 

Reach & Rho (1996) also observed bright [Ol] 63/im emission, particularly enhanced on the edge 
of the remnant where their cut crossed a radio filament, but everywhere higher within the remnant 
than outside. They quote their peak surface brightness as 10^^ erg cm^^ sr^^ (presumably also 
s~^), and estimate the total luminosity of the remnant in this line at 1000 L©. This line is an 
important coolant in shocked gas, but after an unsuccessful attempt to model it as shock emission, 
they interpreted it as evidence that W44 is interacting with a dense molecular cloud, even though 
the region explored does not appear from the IRAS maps to be abutting a dense cloud, being well 



northwest of that explored by Wootten (1977), for example. In Section B.5, we show that the [Ol] 
63/um emission is in fact quantitatively consistent with our model, involving no interaction with 
dense molecular material. 



As mentioned previously, the OH masers observed by Claussen et al. (1997), were also inter- 
preted as being collisionally excited by the SNR shock passing through a molecular cloud. They 
point out that DeNoyer (1983) contradicted the early evidence for such interaction, but then cite 
the [Ol] observations just discounted as the strongest independent evidence for it. 

Claussen et al. (1997) were further influenced by the presence of an unidentified EGRET 
gamma ray source error circle overlapping the eastern limb of W44, suggesting that it may be due 
to cosmic rays from W44 interacting with an abutting dense cloud. The Wooten cloud is in fact 
projected on that location and Esposito et al. (1996) show that the source brightness, 50 x 10^* 
photons cm~^ s~^ above 100 MeV, would require a cosmic ray intensity within the cloud nearly 
100 times ambient, even if the entire cloud is being irradiated. 

3 Analytic Modeling 

3.1 Overview of the The Basic Model 

In combination, the radio continuum and 21 cm observations provide a persuasive qualitative 
picture of this object. The strong filamentation of the radio continuum is characteristic of that from 
the very thin dense shell on the surface of a radiative remnant. The radio continuum maps, with 
their concentration toward (and sharp edge on) the east, and elongation southeast to northwest, 
give the impression of an emitting surface that is incomplete, somewhat like half an egg shell with 
the open end roughly toward the west or southwest. The 21 cm emission is also suggestive of an 
incomplete shell, with the observable part in recession. 

For our model, we adopt this qualitative picture, assuming that evolution is occurring in a 
medium that is denser on the northeast far side, that shell formation has occurred, but only on 
the denser end, and that the filamentation is indicative of irregularity (corrugation) in the shell's 
surface, as explored for optical filamentation by Hester (1987). A schematic of the model is shown 
in Figure |^. 

Our choice of northeast to southwest, orthogonal to the observed major axis, for the projected 
axis of the density gradient requires a bit of explanation. We had initially expected that the density 
gradient would provide the elongation of the remnant, and that we should be taking the gradi- 
ent along the elongation. When we became more familiar with the Kompaneets Approximation 
(Maciejewski and Cox, 1998), however, we discovered that in an exponential density distribution a 
remnant remains nearly spherical much longer than we had expected. It expands more rapidly in 
directions of lower density, but the effect is to displace the center of the sphere down the gradient 
from the explosion site, not to elongate the envelope. Eventually we realized that a partial shell 
capping one end of a sphere could appear elongated, but that the elongation would be transverse to 
the density gradient. For this reason we have adopted the observed semimajor axis as the remnant 
radius, and taken the model gradient orthogonal to it. 

The remnant characteristics, summarized in Table || can be combined in a variety of ways to 
demonstrate that W44 must indeed be a supernova remnant, with an explosion energy of about 
]^q5i g]-g8_ por brevity we will adopt this result and from it derive the remnant characteristics 
assuming the age to be 20,000 years, the nominal age of the pulsar. There are two free parameters, 
the density at the explosion site and the scale height of the density distribution. 

Determination of these two parameters follows from three requirements: (1) that shell formation 
has occurred on the dense end and extended about half way around the remnant to the equator, 
that (2) there be sufficient mass already in the half shell to explain the amount seen in 21 cm, but 
(3) with a shell expansion velocity as large as that seen. 

In the remainder of this section, we determine these parameters and make analytic predictions 
for the basic structure (the remnant radius, shell expansion velocity, edge temperature, central 
pressure, central density and temperature, shell density and magnetic field), and for anticipated 
characteristics in observed emission modes (thermal X-rays, radio synchrotron. Ha, 63 /im, and 
gamma rays) . 



V^ildi cLL Ucl lb Lie 


Value 


wubei vauion 




SNR Angular Size 


25 X 35 


Radio Continuum 




Observed Shape 


Pear like 


Radio Continuum 




Elongation 


Northwest-Southeast 






Evolutionary Stage 


Partial Shell Formation 


HI Emission 


1 


Model Axis 


Northeast-Southwest 




2 


Distance 


2.5 to 2.6 kpc 


HI, OH Absorption 


3 


Radius 


~ 11 to 13 pc 






Age 


~ 20, 000 years 


Pulsar 




Radio Continuum Flux Density 


~ 500 Jansky @ 100 MHz 






Spectral Index 


- 0.33 to 0.4 


Radio Continuum 




Shell Expansion 


150 ±15 km sec 


HI Emission 




Shell Nh 


~ 3 X 10^9 cm-2 


HI Emission 




Optical Configuration 


filamentary, within radio remnant 


Ha, [SII] 




jreaK ncx suriace uiigiiLiieos 


~ ± X lu eigb cni s arcsec 






Peak 63 /im surface brightness 


10~^ erg cm^^ sr~-^ 


ISO 




63 iim luminosity 


~ 1000 L© 






Gamma-Ray Source Box 


includes eastern limb 






Gamma-Ray Flux above 100 MeV 


50 X 10~^ photons cm~^ s"^ 


EGRET 




X-Ray Morphology 


Centrally Bright, No Shell 


0.2 to 4 keV X-ray 


4 




Extends Dimly to Radio Boundary 




5 


X-Ray Emission Mechanism 


Thermal 


Observable Lines 




Central Temperature 


- 4 - 8 X 10^ K 


X-ray 




Interior 


~ 0.4 to 0.75 cm-3 


X-ray 


6 


Central Pressure 


1.4 X 10~^ dyn cm~^ 


X-ray 




Absorption Nh 


~ 1.6 - 2.1 X 10^2 cm-2 


X-ray 




Pulsar Nebula 


Unimportant in 0.2 to 2 keV 


4 - 9.5 keV, 20 cm 





Table 1: Quantities derived from the observations: Note 1: The dense shell formation is more 

fully developed on the northeast rear section; Note 2: The model axis is assumed perpendicular 
to the elongation in radio continuum emission because the shell is partial; Note 3: This is a new 
distance estimate. Note 4: Harms et al. 1996 present a 0.5 to 4.0 keV image which shows the 
central brightening, as does Rho et al.^s (1994) ROSAT data. Note 5: The dim envelope of the 
x-ray emission extends to the boundary in most directions. Note 6: The reported electron densities 
have been scaled by 1.07 to reflect a change in the assumed distance from 3000 pc to 2600 pc. 



3.2 Review of Shell Formation in a Homogeneous Medium 

In our notation, n and rie refer to the densities of nuclei and electrons, respectively. Then n = rin + 
riHe = 1-1?^H, where Uh and tih^ are the densities of hydrogen and helium nuclei, and n^^/n^ = 0.1 
is assumed. In this notation, the mass density is p = mn, where the average mass per nucleus is 
m = (1.4/l.l)mH, TOh being the mass of hydrogen. The perfect gas law, assuming equal electron 
and ion temperatures is then p = xnk^T, where X = (^^ + ne)/n. For neutral gas, x is 1, for fully 
ionized it is (2.3/1.1). 

The Sedov self similar structure and evolution yield several useful results for middle aged rem- 
nants if the effects of cooling have not yet been large. The shock radius is Rg = [2.025Eot'^ / (mrio)]^/^ 
where Eo is the explosion energy and Hq is the ambient value of n, assiimcd uniform. The shock 
velocity is then Vs = dRs/dt = (2/5)i?s/i. From the jump conditions, the post shock density is 
Auo, the pressure is ps = [3/4]mnow|, and from the perfect gas law in our notation, the post shock 
temperature is Tg = (3/16)mt;|/(xfcB)- In addition, the pressure in the center of the remnant is 
Pc ~ O.Slps- 

Kahn (1975, 1976) introduced a very useful approximation technique for estimating the cooling 

timescale for supernova remnants. In that approximation, the cooling coefficient takes the form 
L = aKahnT~^/'^ . It has the special property that in the absence of thermal conduction, the time 
for a parcel of hot gas to cool from an initial state is independent of its history and equals 

AUo, = ^. (1) 

For a newly shocked parcel of gas. p, T, and n arc known and thus, so is Afcooi- To find the actual 
time ^finish at which the parcel will finish cooling, the post shock cooling time must be added to 
the time tghock at which the parcel is first shocked (^finish = ^shock + Atcooi)- To find the time at 
which the first parcel cools, the minimum is found for the function ifinish(^shock)- This minimum 
is the cooling time for the remnant, tcooi- Details appear in Cox and Anderson (1982) and Cox 
(1986). 

The cooling time, however, should not be confused with the shell formation time fshcii- The 
first gas which cools in a remnant does so rapidly. It drops to a very low pressure relative to its 
surroundings, and is only later compressed and repressurized by flows of cooling gas from both 
sides, these driven by the pressure differential. It is therefore not correct to assume that "cooling" 
and "compression" (i.e. dense shell formation) proceed simultaneously. From experiments with 
hydrocodes, we estimate that tgheii ~ (7/6)tcooi- The result, and numerical forms with aKahn = 

1.3 X 10-1'^ cgs, are 

teoo, « 0.208 [(m9/"i;3/2)/[^5/2^^(^^^)5/4„4]]Vr 

icool « 45.5 X 10^ years i;5\/'Vny^ (3) 
Ueii « 53 X 10=^ years £;5'("/ny^. (4) 

We are now in a position to estimate the remnant's dynamical characteristics at the moment 
when shell formation reaches the equator, and to make an initial estimate of the ambient density, by 
solving for the epoch at which tshoii(£^o, ^^avg) equals the remnant age, about 20,000 years. Including 
the possibility that the age or explosion energy are probably not exactly our nominal values, the 
values for shock radius, shell velocity, and central pressure when the shell formation time equals 
the age, and the required average density are (with t2o = t/{2Q, 000 yr)) 

i?s « 12pci?5i^/*'t2o'/' (5) 

t;sheu « 0.75fs«176kms-i.B5i'/V*2o'/^ (6) 

Pc w O.Slps ^ 1-5 X 10-9 dyn cm.-^ E^i^^^ /t2o^/^ (7) 

no « 5.5 cm-3i;5i^/V*20^/^ (8) 



These are the (approximate) anticipated conditions at the moment when shell formation first 
occurs, a condition which we believe applies at the equator of W44. The remarkable fact is that 



even at the nominal values of explosion energy and age, the radius and central pressure are easily 
within the observational uncertainties. The shell velocity is somewhat high, but we expect a lower 
velocity toward the denser end where the shell is more completely developed, and which in our 
projection is closer to the line of sight. Further insight derives from evaluating the temperature 
that would have been at the edge of the remnant had there been no cooling, considering it to be 
an estimate of the temperature of hot gas at the outer edge of the hot interior (just interior to the 
dense shell) at this epoch. The result is Touter ~ 7.6 x 10^ K E^^^'^lH^i'^'^- As found by Harrus 
et al. (1997), the remnant edge at this epoch is somewhat too cool to provide a strong shell source 
of X-rays. 

There is another consistency check in the HI column density, about 3 x 10 cm" -2 from Table 0. 
If all the mass in a sphere is swept into a thin shell at the edge, the perpendicular column on one side 
should be UoRs/S. This is about 7 x 10^^ cm~^ for our parameters, a consistent value considering 
the uncertainty in the observations, the recency and incompleteness of shell formation, and the 
unknown fraction of the mass that is molecular. If the observed HI arises from any other cause, it 
would be exceedingly unlikely for it to have the column density and expansion velocity expected 
for the cooled SNR shell. Making the details of the shell mass work out correctly depends only on 
making the correct choice for the density scale height, a topic we defer until after considering the 
effects of thermal conduction on the central density and temperature of the remnant. 



3.3 The Central Density and Temperature 

As we have suggested in the introduction and shown elsewhere (Slavin and Cox, 1992, Cui and Cox, 
1992, and Shelton, Smith, and Cox, 1995), thermal conduction in the hot interior of a remnant 
prevents formation of the vacant interior characteristic of nonconductive evolution. Transporting 
heat away from the remnant interior lowers the specific entropy and therefore the expansion rate 
needed to accommodate the falling pressure, thereby moderating the extraordinary reduction in 
central density found in models ignoring conduction. The central pressure remains uniform at 
about 0.3 of that at the edge, but the temperature and density are also made nearly uniform as 
well. In this section we show that this effect can provide the conditions found in the x-ray emitting 
interior of W44. 

There is considerable concern, of course, that thermal conduction might be very heavily sup- 
pressed, and/or anisotropic in supernova remnants due to the magnetic field. We have addressed 
this point previously (Slavin and Cox, 1992, and Cox, 1996) but, in light of its importance to our 
present application, discuss it somewhat further in section 



4.4 



Experience with hydrocodes has suggested that the central density Uc at the epoch of shell 
formation, in a remnant with unsuppressed conduction, is roughly O.Itt-q. We have previously shown 
(Shelton, Smith, and Cox, 1995) that this factor of 10 density contrast between the interior and 
the ambient medium can be understood analytically, and elaborate on that result now, exploring 
its dependence on the magnitude of the thermal conductivity. 

Cui and Cox (1992) developed a similarity solution for conduction that made it possible to 
estimate a remnant's central adiabat, temperature, and density analytically for times in the vicinity 
of the cooling epoch. There was unfortunately a typographical error; the density exponent in their 
equations (A17) and (39) should be 208/225. In addition, the units are incorrectly described just 
before (39) but are correct in their appendix. The generalization of the coefficient B in their 
equation (A14), to find the dependence on the magnitude of the conductivity, is 

B = 29.2/3ml2/25„208/225/[(^^^)7/2^26/75] (g) 

where the classical conductivity is taken as /ST^/^, with f3 ^ 6 x 10^^ cgs. 

With this result, and assuming the Kahn cooling function, the ratio of the central to ambient 
density turns out to be independent of explosion energy or ambient density. The value and form 
are, 

njn, = 0.334gi8/7i(i_j/^)i8/25 [(t/i_,)i/25 _ (o.44g^/i4)i/25] (^q) 

with q = lAdmaKahnP / ixkn)^ (defined so that q = 1 for full conduction and normal cooling rates). 
For t — (7/6)tcooi, as at shell formation, and q = 1, we find nc/rio « 0.13, nicely in agreement 



with hydrocode experience. Equation plj| probably does not successfully describe time dependence, 
but should provide a guide to the dependence on conductivity. It predicts the central density to 
be proportional to (3'^-'^^^. A factor of 10 reduction in conductivity should lead to only a factor of 
1.8 reduction in the central density and increase in central temperature. 

With rio ~ 5.5 cm^'^, and q ^ 1, we expect a central density within the remnant cavity of 
ric ~ 0.7 cm"'^, just the level found from the x-ray surface brightness of W44 (compare with the 
x-ray results for the central density in Table |l|) . As the central pressure was shown in the previous 
section to be at the correct level, the temperature will be appropriate also, roughly 7 x 10® K. 
We can therefore be confident that conduction will lead to x-ray emission similar in characteristic 
temperature and emission measure to that seen. We defer to the hydrocodes, however, to learn 
whether the spectral details in the presence of nonequilibrium ionization, and the radial profile of 
the surface brightness are also well approximated. 



3.4 The Radio Continuum 

3.4.1 Evaluation of The Contribution of the van der Laan Mechanism 

There is considerable interest in the possibility that the bulk of the cosmic rays arise from diffusive 
shock acceleration in supernova remnants, an interest which is encouraged by their powerful non- 
thermal emission (Ellison et al. 1994). Before assuming that diffusive shock acceleration plays the 
dominant role in an evolved supernova remnant which is bounded by a radiative shock, however, it 
is important to consider the synchrotron emission that would arise from swept up ambient cosmic 
ray electrons, adiabatically compressioned into the dense shell and radiating in the compressed 
ambient magnetic field (van der Laan, 1962). 

Neglecting small inefficiencies due to misalignment and isotropization, a co-moving relativistic 
particle in a region compressed by a factor x would have its energy betatron-enhanced (2d adiabatic 
compression) by about a factor of x^^'^ . The net effect for a supernova remnant with a cooled shell 
of mass M and original density po is that cosmic ray electrons previously occupying a volume 
Mj Pa of interstellar space find themselves compressed into a shell whose volume is x times smaller, 
enhanced in their energies by a factor of x^^^, and radiating in a magnetic field about x times 
higher. The characteristic frequency of synchrotron emission of an electron of energy i? in a 
field of strength B is proportional to and is therefore a factor of x^ higher. Similarly, the 

total power per electron, proportional to (BE)^, is raised by a factor x^. For a distribution of 
electrons radiating synchrotron with an emission spectrum Av^" per electron before compression, 
the spectrum per electron after compression is with A' = a;^+^"A, for a total gain at fixed 

frequency of G{x) = a;^+^". 

If one is willing to make assumptions about the isotropy of the particle distribution before, 
during, and after the acceleration, then this estimate can be made somewhat more precise. For 
example, assuming an initially random distribution of pitch angles and the post acceleration re- 
randomization, the overall gain per electron from van der Laan (1962) is 



G{x) = X 



l+a 



2^2{x~iy/^ ^ X > 



(11) 



Compared with the previous estimate, the numerical difference is slight (here G{x) k, (tt / A)"^" x^'^'^" 
at large x). 

Therefore, except for the required units conversions, the flux density of the remnant's radio 
continuum from this mechanism is 

S{,y) = G{x)Voea{i^)/d^, (12) 

where Vo — Mghcii/Po is the shell volume before compression, eo{i>) is the synchrotron emissivity 
of the ambient ISM at frequency x is the shell compression factor, and G{x) is the gain as given 
above. The next few paragraphs assemble these parameters. 

^From the observations summarized in Table |l| and the subsequent modeling, we have estimated 
the shell radius at 11 to 13 pc, the shell formation process complete over half that surface, the 
one sided surface density at 3 x 10^^ cm~^, and the initial density at 6 cm~'^. Thus, the mass of 
material now in the half shell is about 300 Mq and its initial volume about Vo — 1600 pc"^. 



W44 has been studied between about 10^ and 10'* MHz. Unfortunately, prior to their being 
swept up by the remnant, the responsible electrons radiate only below 10 Mhz where their study 
is difficult. Thus, the ambient emissivity must be evaluated at the lowest observable frequencies 
and then extrapolated into the frequency range observed in the remnant. Webber et al. (1980) 
and Rockstroh and Webber (1978) conclude that the data available appear to be consistent at the 
lowest frequencies with an ambient interstellar synchrotron emissivity of roughly 

eo - 1.5 X IQ-^" [i//30(MHz)]-^^ W(m^sr Hz)-\ (13) 

with apparently real place to place variability. Extrapolation to 100 MHz then yields 

eo(100 Mhz) - 0.8 X 10"''° W(m^sr Hz)"^ (14) 

The fact that the emission spectrum quoted above, with a ~ 0.57, is steeper than the observed 
spectrum for W44 requires further elaboration. It is attached to the conclusion of Webber et al. 
(1980) that the ambient electron particle index {l+2a) is 2.14 between 70 and 2500 MeV. But 
they further conclude that below 70 MeV the particle spectrum flattens toward an index of 1.6. In 
the shell of W44, the emission at 100 MHz derives from electrons with energies about 250 MeV, 
but prior to compression into the shell the energies of the same electrons would have been 7 to 
10 times lower, or roughly 30 MeV. As these electrons are in fact in the flattened part of the 
spectrum, a spectral index lower than 0.57 but greater than 0.3 (corresponding to particle index 
1.6) is expected. Therefore the observed spectral index of W44, 0.33 to 0.4, is nicely consistent. 
There should be a gradual steepening of the radio spectrum at high energies, but the observability 
is complicated by the distribution of compression factors within the remnant shell. 

This spectral curvature complicates the simple analysis we have undertaken. We have elected 
to use an intermediate effective index, acff — 0.5, in evaluating the gain, both because it represents 
the overall increase in the number of suitable radiating electrons better than the observed indices 
of W44 or the ambient ISM, and because it leads to a very nice summary formula for the predicted 
flux density. We expect this simplification to provide an estimate which is reliable to within a 
factor of two except for remnants with very flat radio spectra. Our adopted approximation for the 
gain is therefore 

G{x) « Jx^ (15) 

In attempting to estimate the compression factor in the shell, we make use of the fact that 
with either analytic approximations or hydrocode models, it is the total pressure in the dense shell 
that is its best known characteristic. Then, making assumptions about the thermal, magnetic, and 
cosmic ray contributions, the corresponding compression factor can be determined. 

In the Sedov solution, the pressure at the densest part of the edge is about 3.3 times higher than 
the central pressure. When cooling sets in, the thermal pressure near the edge drops precipitously 
at first and then recovers somewhat as neighboring material compresses the cooled gas. On the 
dense side of W44, the pressure in the dense shell should have recovered in this way, but to a 
somewhat lower value than for a Sedov remnant which had had no cooling. Cox (1972) suggested 
that the shell velocity would probably be about the same as the post shock velocity just before 
cooling, leading to a shell pressure that is 3/4ths of the precooling edge pressure. With this 
guidance, we estimate that the shell pressure in W44 is now about twice the central pressure pc, 
where from Table |l|, we have pc lA x 10~^ dyn cm~^. 

For an overall compression factor a;, the thermal pressure in the dense shell will be (xks'>T'oTf)x, 
where Tf is the temperature in the shell. The final temperature is not well known, but assuming 
Tf ^ 1000 K, and x = 1 provides a reasonable estimate in the neutral gas. With rto ~ 6 cm"'^, the 
thermal pressure in the shell is negligible compared to the nonthermal terms treated next. 

In the shock and subsequent compression into the shell, it is expected that the tangential com- 
ponent of the magnetic field will increase in proportion to the density. As a rough approximation 
(also used in our hydrocodes), we take the field as effectively transverse everywhere, in which case 
the magnetic pressure is x'^B'^_^/{87r). We consider values of Bo± in the range 2 to 5 /iG to be 
plausible. 

The pressure due to the compressed ambient cosmic ray electrons is negligible, but the ions could 
well be important. Their ambient pressure is comparable to that of the magnetic field and their 



non-relativistic component will have its pressure increasing approximately as , like the magnetic 
pressure, so long as they do not isotropize during the compression. Thus, their contribution to the 
shell pressure, after strong compression, could be comparable to the magnetic pressure. It can be 
substantially less if the cosmic ray ions are not efficiently localized (i.e. they leak out) or greater if 
they are accelerated by more than just the compression. In what follows, the cosmic ray pressure 
is taken as pcR = /?Pb, where /3 is expected to be in the range to 1. 

The total pressure is thus approximately (1 + j3)x'^ B'^^/ {^i:) which must equal Pshoii, yielding 
an estimated compression factor of 



90 



(l+/3)0-' 



Br, 1 



(16) 



which ranges from 40 to 130 over the ranges of (3 and Bq. The implied range of shell density is 220 
to 800 cm~^. The magnetic field within the shell is less uncertain, ranging from 185 to 260 /iG, 
depending only on (3. 

This result can be generalized to any remnant with known shell expansion speed and preshock 
density, realizing that what we have done is to equate the nonthermal pressure in the shell with 
the ram pressure of the radiative shock, Pram = Po'f^shcii' yielding 
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Thus, with Q!cff — 0.5, the gain can be written 
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the initial volume of material in the shell as 
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and the resulting synchrotron flux density as 
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For comparison, Table |l| shows the measured flux density for W44 at 100 MHz to be about 500 
Janskys, within the uncertainties, in complete agreement with this estimate. 



3.4.2 Discussion 

Our estimate of the magnetic field in the shell of Sghcii ~ 200 /iG was made before learning of 
a measurement of the field in W44 by Claussen, et al. (1997), who also found of order 200 /iG. 
Their measurement was of polarization of OH masers, was more uncertain than our estimation 
procedure, and was attributed to quite a different picture of the source region. It is therefore 
somewhat less satisfying than the numerical coincidence would indicate, but could be appropriate. 
Their interpretation is that a very slow shock is being driven into an abutting molecular cloud, 
with a density close to lO'* cm~'^ and an ambient field not much smaller than what they measure. 
But the close association of the OH masers with apparently tangentially moving regions of radio 



continuum emission and the fact that the measured fields can be achieved in the remnant shell 
cause one to question whether a molecular cloud interaction is necessary. Certainly the small 
compressions involved in their scenario would yield negligible radio synchrotron emission via the 
mechanism we are investigating. 

For densities in the range found for the shell, 220 to 800 cm^^, the recombination time is only a 
few hundred years, so it is reasonable to expect the dense shell to be mostly neutral and responsible 
for the observed 21 cm radiation. 

With a shell compression factor of roughly 100, the local shell thickness (before any distortions) 
is less than half a percent of the radius. This implies that the radio synchrotron emission is very 
highly edge brightened and subject to the high degree of ripple/tangency filamentation described 
by Hester (1987). We present some simulations of this effect in section 5.2.3 . 

The synchrotron emissivity estimate above hides a fundamental problem with this approach. 
We have assumed that the best estimates of the ambient cosmic ray electron population, the syn- 
chrotron emissivity, and the tangential field strength are mutually consistent. As discussed by 
Rockstroh and Webber (1978), the characteristic field strength required for consistency is surely 
greater than the mean field, with the rms field being more suitable. Results of attempts at con- 
sistency depend somewhat on the assumed spectral index, but typically require a characteristic 
field of about 9 /iG. The value of Eq appropriate for a 3 /iG field should then be reduced (for 
tteff = 0.5) by about a factor of 5. Therefore, the flux density of W44 cannot be matched by the 
van der Laan mechanism alone unless the ambient CR electron population is about a factor of 
3 to 5 larger, depending on the shell mass, than estimated for the solar vicinity. (Alternatively, 
the anisotropy of the synchrotron radiation formulae may have fooled us). We will return to this 
subject in discussing the gamma rays anticipated from the shell. 

The radio continuum flux predicted from the van der Laan mechanism is similar to that ob- 
served, but there is clearly room for contributions or enhancements by additional processes, so 
long as those processes do not produce a markedly different spectrum. The next three paragraphs 
discuss some of the possibilities. 

Some injection by the shock of new cosmic ray electrons may be occuring. Alternatively, 
electrons accelerated by the remnant when it was younger, or by the pulsar, may be caught up 
in shell formation and forced to radiate in the vastly enhanced magnetic field there. One might 
suppose that the problem with either possibility is that there is no obvious reason for the surface 
brightness and spectrum to be so close to those found for the van der Laan mechanism and by 
observation, suggesting that we consider, as below, factors that could enhance the radio synchrotron 
emission over the above estimate without invoking a totally independent mechanism. But perhaps 
there is a reason for achieving approximately this result with other mechanisms, namely that our 
emissivity is basically the cquipartition result-if for a given pressure, ambient field, and shell mass 
we ask for the emission from the usual percentage of electrons associated with cosmic ray ions 
in cquipartition with the shell's magnetic field, we will have roughly what we have calculated, 
independent of the details of the acceleration mechanism. 

Ambient cosmic ray electrons may receive substantial diffusive acceleration on passing through 
the adiabatic gas subshock, followed by betatron acceleration during compression into the dense 
shell. These ambient electrons, apparently having a slightly shallower index than that of newly 
injected and accelerated ones, will maintain their original spectral index while achieving an average 
energy boost per particle. From elementary theory, the magnitude of the boost is unknown, but 
potentially large. 

Shell formation offers yet another acceleration possibility. Electrons which attempt to diffuse 
away from the cold dense shell find themselves continually swept back in because there is ap- 
proaching flow on both sides: recently swept up ambient material cools rapidly and is added to the 
shell from the outside, while material just interior to the shell which was heated long ago finally 
cools and provides a condensation fiow to the inside (occasionally accompanied by another weak 
radiative shock) . This two sided fiow convergence not only reduces particle losses from the densest 
region, but offers an as-yet-unexplored site of first order Fermi acceleration. 



3.5 The Ha, [Ol] 63/jm, and Related Luminosities 

A rough estimate of the Ha luminosity of the remnant can be made following the methods of Cox 
(1970, 1972a). There are several factors to consider. Hot material inside the remnant and newly 
shock-hcatcd material on the outside edge arc cooling off, rccombining, and joining the dense shell. 
They contribute directly to Ha, but also emit EUV which fluoresces in the neighboring neutral 
material, shell and ambient, to yield even more Ha. A further complication is that fluorescence in 
the ambient medium is not instantaneous; ionization is effectively stored for a recombination time, 
which at rio ~ 6 cm""^ is about the age of W44. 

Because the recombinations in shell and ambient material can be distinguished via both spatial 
distribution and in velocity, we consider them separately. Following Cox (1972a), we estimate 
that in cooling from high temperature, approximately 4 ionizing photons are radiated per atom. 
Assuming that half of these photons escape to ionize the ambient medium, the ionization rate for 
the surroimdings is 2dN/dt and the recombination rate for the shell is idN/dt (there being one 
extra from the recombination of the cooling atom). (Here dN/dt is the rate at which atoms join 
the shell.) The total number of atoms in the remnant is approximately UoV 10^° (about 1000 
M0 ). Of these, we expect that roughly a third (but see below) have cooled during the ^ 3000 
years between tcooi and fgheib contributing 4 x 10^* atoms s~^ to dN/dt. The rate at which atoms 
are overtaken by the radiative shock (presently covering half the remnant's surface, of total area A) 
is na{A/2)vs ^ 1 X 10"*^ atoms We thus expect a total of dN/dt ^ 5 x 10^* atoms s~^. Results 
from our 2d hydrocode (below) show this estimate to be about a factor of 2 too high at present 
(20,000 yr), with the higher rate more representative of the recent past. Because the ionization 
level of the ambient medium averages the rate over a siibstantial time, we take the higher rate for 
it, yielding 1 x 10^^ photoionizations s^^ for the last 3000 years. For the shell fluorescence the 
instantaneous value of dN/dt is more appropriate, yielding 0.8 x 10*® recombinations s~^ in the 
dense shell. 

/,From recombination theory there are approximately 0.6 eV of Ha emission per average (Case 
B) recombination. The shell luminosity in Ha from the radiative shock should thus be about 
3 X 10^^ erg s"-*^. From the above analysis, however, the total number of atoms presently joining 
the dense shell, including those from the cooling interior, is 2 to 3 times higher, implying a total 
Ha luminosity of about 7 x 10^^ erg s~^. Its spatial distribution and velocity structure should 
resemble those of the radio continuum and HI shell respectively, with one caveat: because shell 
formation is presently occurring at the equator of W44, we might expect the emission to be brightest 
there, rather than in the east or northeast where shell formation has already passed its peak. This 
equatorial configuration is in fact apparent in Fig lb of the Ha picture of Giacani et al. (1997) (who 
attributed it to enhanced absorption in the east), while their Figure 5 demonstrates the expected 
strong relationship between the Ha and radio continuum enhancements and filamentation. 

Our luminosity estimate corresponds to an average surface brightness over the ^ 500 square 
arcmin of the remnant of 6 x 10~^^ ergs (cm^ s^ arcmin^)^^. The extinction at Ha implied by an 
x-ray absorption column density of about (1.6 to 2.1) x 10^^, however, is a factor of 600 to 4400, 
for an observable average of (0.1 to l)xlO~"'^^ ergs (cm^ S"'^ arcmin^)^^. Giacani et al. (1997) 
quote their peak observed surface brightness as 10~^^ ergs (cm^ s^ arcmin^)^^, with the dimmest 
observable features at 10"^'' ergs (cm^ s^ arcmin^)"^. Their brightest features, a factor of 10 to 
100 higher than our average, are strongly c^nhanced by filamentation; but their dimmer features, a 
factor of 1 to 10 above our estimate, appear much smoother in the equatorial belt. Given the fact 
that we expect concentration in that region, and that the large extinction is not only uncertain 
but likely to be variable, the agreement is as good as can be expected. 

Because the recombination time in the ambient medium is long, its luminosity will be lower. 
There have been a total of roughly 10^'" ionizing photons captured in the surroundings. Neglecting 
recombinations, they will ionize a similar number of atoms, which will then recombine with a 
timescale of 10"'' years /no ~ 1.8 x 10'' years, yielding 0.2 x 10''^ recombinations s~^ and a luminosity 
of 2 x 10^^ erg s^^ . Although the total radiation by the shell is expected to be about 3 times brighter, 
and enormously higher in surface brightness in restricted regions, the fluorescence in the ambient 
medium is not inconsequential. The ionized mass should be about 1000 , like the mass within 
the remnant, and extend roughly 3 pc outside (even farther but dimmer on the low density end), 
with an emission measure exceeding 100 cm~^ pc. In the absence of inhomogeneity other than 



the density gradient, it should be brightest behind the northeast part of the remnant, where the 
ambient density is high. 

Of course the emission of the cooling gas is not restricted to Ha. The overall spectrum should 
be rather similar to that of a radiative shock with a speed of about 150 km s"'^. From Model 
G of Raymond (1976), for example, we expect emission from optical lines of [OIII], [Oil], [Ol], 
[Nil], and [SII], and near IR lines of [Cl] that are comparable in intensity to Ha. In the infrared, 
predicted strong lines are [Fell] at 26 /xm, [Ol] at 63 /xm, and [Nell] at 12.8 /xm, the first of which 
is about twice as bright as Ha, the others about one third. The basic picture is not much different 
in more recent calculations (although inclusion of proton impact excitation alters the strengths of 
the infrared lines somewhat). Essentially all shock models in the velocity range 130 to 160 km 
sec^^ yield 1 to 3 eV of Ha and 0.5 to 1 eV of 63 fim emission per average atom through the 
shock. Once again assuming that the radiation by material joining the dense shell from inside is 
similar in spectrum to the radiative shock, we conclude that the luminosity in any of the above 
lines should be about (3 to 10)xlO'^^ erg s~^. The [SII] was observed by Giacani et al. (1997), 
with a distribution (and presumably intensity, though it was not discussed) similar to that of Ha. 
The 63 /im line was measured with the ISO LWS by Reach and Rho (1996); from the regions 
they sampled they estimate a total luminosity of about lO'^L©, in excellent agreement with our 
estimate. 

Interestingly, Reach and Rho considered a radiative shock model, but rejected it as unable to 
provide the observed emission. There appear to be three contributing reasons for our difference: 
they compared the theory of a single shock viewed face on with their maximum observed surface 
brightness, 10~^ ergs cm^^ s~^ sr~^, rather than an estimated average which for IQ^Lq would be 
a factor of 10 lower; in addition, the average surface brightness of a spherical emitting surface is 
four times that of the local normal value for one side; and, finally, we have an additional factor 
of 2.5 to include the effects of current shell formation. The net effect is a factor of 100 which is 
just the difference between their estimate of the rate per unit area that the shock would have to 
be encountering ions and the 10* atoms cm~^ s~^ in our model. This caused them to believe that 
an unacceptably high ram pressure was somehow present where the [Ol] emission occurred, and 
that that was somehow connected to an encounter with dense molecular material. In our model 
this is unnecessary: the correct total emission is found, and like the radio continuum and optical 
emission, there are local enhancements over the average of more than a factor of 10 due to rippling 
and edge brightening of the thin shell. 



3.6 The Gamma Rays 

It is fairly straightforward to estimate the maximum gamma ray intensity from cosmic ray ions 
interacting with the mass in the dense shell. Compression of the nonrelativistic component of 
ambient CR ions, or acceleration of ions in the remnant can potentially lead to a CR pressure in 
the shell comparable to the magnetic pressure, their sum being the shell pressure, approximately 
2.8 X 10~^ dyn cm^^. This amounts to an enhancement of the ion pressure by a factor of 1400 
over ambient. From the review by Bloemen (1989), and models in Esposito et al. (1996), we infer 
that this enhancement would lead to a gamma ray production rate above 100 MeV, from 7r° decay, 
of about 2.2 x 10^^^ photons s^^ per H atom in the shell. For a shell mass of about 450 M© , as 
found in our 2d modeling later in the paper, and a distance of 2.5 kpc, the maximum flux would 
be 12 X 10^* photons cm~^ s~^, about a fourth of the observed value, concentrated in the east or 
northeast. (The total mass interior to W44, in our model, is about 1000 M© , but in our simplest 
picture we do not expect the more diffuse portions in the interior to be subjected to this high 
cosmic ray pressure. With some leakage and the associated additional compression, however, more 
of the remnant's mass may be exposed.) One may suppose that cosmic rays are leaking out of the 
remnant altogether, and irradiating the nearby dense molecular material, though if they do so at 
greatly reduced pressure, vastly more mass will be required for interaction. (Esposito et al., 1996, 
found that irradiation of the entire Wootcn cloud would require an average cosmic ray enhancement 
over ambient of a factor of about 80.) 

If diffusion prevents the CR ion pressure from limiting the compression to the low value implied 
above (about 40), the bremsstrahlung component of the gamma rays can be enhanced by an 



similarly large factor. Generalizing on the methods of Section 3.4.1, the ratio of Gamma ray 



bremsstrahlung to radio synchrotron enhancement factors turns out to be x"""*'', where x is the 
compression factor and cteff is the effective spectral index, which we have taken as 0.5. Comparison 
with our previous results for the van der Laan mechanism yields an overall gain of 

G^{x) w 800(a;/100)^ (24) 

Again following the review by Bloemen (1989), and models in Esposito et al. (1996), we infer that 

this enhancement, for x = 100, would lead to a gamma ray production rate above 100 McV of about 
6 X 10~^^ photons s~^ per H atom in the shell. For the shell mass and distance given above, the flux 
would be about 4 x 10~^ photons cm~^ s~^, about 30% of the upper limit from -k° decay and 7% 
of the observed gamma ray flux. Recalling from our discussion of the synchrotron radiation that 
the observed flux required approximately 4 times more CR electrons than anticipated from solar 
neighborhood estimates, and noting that electrons in the same energy range would be responsible 
for both the synchrotron and bremsstrahlung, we conclude that the gamma ray estimate must be 
raised by a similar factor, to 16 xl0~^ photons cm~^ s~^, comparable to our upper limit from 7r° 
decay. Both mechanisms therefore require more target mass than our model provides, by about 
the same factor. 

An alternative approach employed by deJager and Mastichiadis (1997), solves for the CR elec- 
tron populations needed to produce both the synchrotron and bremsstrahlung. Taking the ratio, 
they obtain a relationship between the required magnetic field and density in the interaction re- 
gion if both are to be produced. Interpreting the interaction region as the shell in our model, 
their result for the CR population required to produce the synchrotron is consistent with our 
previous conclusions. Their result for the required shell density for generation of the gamma 
rays can be summarized as ridjM « 25 cm~^(i?/10/zG)^-^'^, while in our model we anticipate 
Tisheii ~ 20 cm~^(S/10/uG). At the shell fields of 180 to 260 /xG of our model, the gamma ray pro- 
duction would require a shell density about 7 times that available. We repeated their analysis with 
a slightly steeper index and found the disparity to be slightly less, but still a factor of 4, consistent 
with our conclusion above that the electrons required in our shell to produce the synchrotron fail 
by about a factor of 4 to produce the observed gamma rays. 

The contribution expected from inverse Compton radiation was also calculated by deJager and 
Mastichiadis, assuming a radiation field that is enhanced by the glowing dust in the neighboring 
molecular cloud. They concluded that the contribution was less than that from bremsstrahlung 
when the mean density in the interaction region exceeds about 10 cm~"^, a very low value compared 
to that found in our shell. Inverse Compton will also arise from CR electrons within the remnant 
but not yet swept into the dense shell, but that too appears negligible. 

The net result of this investigation is that the dense shell in our model should produce gamma 
rays at roughly 25 to 50% of the rate actually observed in the vicinity by EGRET. 

3.7 Evolution in a Density Gradient 

Approximation to Kompaneets 

In planning the use of a 2d hydrocode to provide a model including the partial shell formation 
of W44, we assumed the ambient density to have the exponential distribution: 

Po(z) = /9*e-^/\ (25) 

where is at the explosion site, and z is the distance from it, measured down the density gradient. 
Thus, to specify a model, the ID input parameters (the explosion energy and the ambient 
density po) have to be supplemented by the stratification scale-height h. The parameter space 
for initialization of a run is therefore three-dimensional, with age a fourth parameter. As high- 
resolution simulations in two dimensions are expensive, an extensive numerical survey of that space 
is not possible. Therefore, we searched for appropriate input parameters with the help of analj^tical 
calculations based on the Kompaneets Approximation (Kompaneets, 1960). 

The basic idea of the Kompaneets Approximation is to assume that the post-shock pressure ps 
is uniform over the surface of the remnant, while the evolution occurs in an exponentially stratified 
medium. For a strong non-radiative shock in a 7 = | gas, ps = ^p^vl, where Vs is the shock 
velocity and po is the ambient gas density. With uniform post-shock pressure, the local shock 



speed is inversely proportional to the square root of the local ambient density. Thus the evolution 
of the shape of the remnant can be calculated without first knowing the time evolution. Once the 
shape is known, as a function of some size variable such as the semi-major axis, the volume can 
be calculated as well. Then, because the pressure is inversely proportional to the volume, the time 
can be found from an integration. 

Maciejewski and Cox (1998) extend this approach in two ways, the first of which is to assume 
that the structure of the remnant is closely approximated by a prolate ellipsoid of semi major and 
minor axes a and b. They then readily find the semi-minor axis b, and the distance s by which the 
center of the ellipsoid has displaced from the explosion site, as functions of semi-major axis a. A 
surprising result of this exercise is that even remnants extending over several scale heights remain 
nearly spherical; the principal geometrical effect of the gradient is to shift the remnant center 
gradually away from the explosion site. (See also Dohm-Palmer and Jones, 1996, for a thorough 
discussion of this behavior based on 2d hydrodynamic modeling.) 

Assuming the postshock pressure, explosion energy, and volume have the same relationship as 
in the ID Sedov solution determines the dependence of the shell expansion velocity on the SNR 
volume, and an expression for the time derivative of the semi-major axis. Numerical integration 
then yields the time evolution t{a), at which point, all of the quantities parameterized by the 
semi-major axis a can be expressed as functions of time. 
Cooling and shell formation times 

Maciejewski and Cox (1998) also evaluate the progression of the cooling time over a remnant 
in a density gradient, in almost exact parallel to the 1-d case. With the Kahn approximation, the 
incremental cooling time Aicooi of a recently shocked parcel depends on the post shock pressure 
and the local preshock density. Adding that to the time at which a parcel is shocked determines 
its time to finish cooling. Minimizing tfinish over parcels (for a particular direction in the 2d case) 
determines the cooling time on that part of the remnant. Finally, multiplying by 7/6ths supplies 
an estimate of the time for shell formation. By performing this calculation for the remnant's dense 
end, tenuous end, and equator, the combination of Kahn and Kompaneets approximations provides 
a useful tool for the exploration of shell formation over the surface. 

Even a modest density gradient leads to a great spread in cooling and shell formation times for 
different parts of the remnant's surface, making it possible to select input parameters which will 
lead to any desired degree of partial shell formation. 

We examined a wide range of input parameters. In particular, we used our analytical method 
to generate plots on the p^-h plane for assumed and linear size of the remnant. The plotted 
quantities included the time required to reach the assumed size, the post-shock pressure, the cooling 
and shell formation time-scales for the dense, equatorial, and tenuous directions, and the shock 
velocities in those directions. (Although the Kompaneets model, like the Sedov model, cannot be 
used very reliably to estimate quantities after shell formation has occurred, we make the usual 
approximation that shortly after shell formation the shell velocity is about 3/4ths of the shock 
velocity predicted for the non-radiative evolution.) In the limit of large scale heights, the results 
are consistent with our ID calculations. 
W44 's parameters 

In this way we were able to locate regions of parameter space which promised to satisfy the 
observational constraints discussed for the 1-d model, adding the further condition that shell for- 
mation must be well underway on the dense end, just commencing at the equator, and not yet 
present on the open more tenuous end. This must be reconciled with a substantial radial velocity 
for the receding part of the 21 cm shell. 

These constraints turned out not to be sufficient to determine the scale height uniquely. In 
the end we chose to examine the most asymmetric (smallest scale-height) case that appeared to 
be viable; in it the cooling times on the high and low density ends differ considerably, but the 
expansion velocity into the densest medium is still large (models with still smaller scale-heights 
have too low an expansion velocity on the dense end). The selected input parameters and estimated 
characteristics of the remnant at 20,000 years are shown in Table |2| while Figure ^ summarizes the 
configuration and goal of the evolution. 



Table 2: 2D Model Parameters and Predicted Characteristics 
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3.8 Sgro-Chevalier and Approximate Scalings 

As it is unhkely that the hydrodynamic models will produce a remnant of precisely the size desired, 
it is useful to be aware of the Sgro-Chevalier scahng law (Sgro, 1972, Chevaher, 1974), that each 
evolution represents a family of evolutions with the same value of the parameter En^. This scaling 
was based on the fact that both ionization and cooling have temperature dependent two-body 
rates, and extends itself to include models with differing electron and ion temperatures, thermal 
conduction, and dust sputtering. In our case, for example, the hydrodynamic runs delivered models 
with radii of 11 pc, which would need to be placed at a distance of 2200 pc to have the correct 
angular size. The Sgro-Chevalier family member with a radius of 12.5 pc and distance of 2500 
pc has the same temperature and velocity structure and the same distribution of ionization, but 
has an initial density of 5.5 cm~'^, an energy of E51 = 1.3, an age of 22,700 years, x-ray surface 
brightness and count rate (at the same angular size) reduced by a factor of 0.88, and unchanged 
synchrotron flux density. The shell mass is higher by the same factor as the energy, 1.3. 

A more general but more approximate scaling, searching for a remnant at the same evolutionary 
stage as a specific hydrodynamic model but with a different radius, can be based on Equations ^, |2^, 
and noticing that for coUisional equilibrium plasmas with temperatures in the range 3 to 8 x 10^ K, 
the ROSAT counts per emission measure function is approximately proportional to T^, particularly 
with significant absorption. The latter implies that the x-ray band luminosity is approximately 
proportional to R^p^ . The inferred scaling is that at a given evolutionary state, the various model 
parameters scale with rio and Enl as: -Rshock oc n~^{En1Y/'^ , t oc n~^{En1Y/^^, v oc (E'ri^)^/^'*, 
T cx {Enlf/\ 5,_Ray (X n,{EnlY/\ F^^n.y « O^noiEnl)^/' , and i^,y„eh cx 0^{Enl)^/', where S 
and F refer to surface brightness and flux, and 9 to angular size. Thus, scaling a model from one 
with a radius of 11 pc to one with a radius of 12.5 pc and constant E^i requires an initial density of 
4.6 cm~'^, and describes a remnant with an age of 23,700 years, x-ray surface brightness and count 
rate (at the same angular size) reduced by a factor of 0.53, and a synchrotron flux density lower 
by 0.77. The shell mass is higher by 1.09, and unlike the Sgro-Chevalier scaling, the velocities 
and temperatures are reduced by factors of 0.96 and 0.92 respectively. With the slightly lower 
temperature the ionization and therefore the details of the spectrum will be slightly different as 
well. 

We will return to these results in discussing the hydrodynamic models. The main feature, 
however, is that all changes are moderate, essentially within uncertainties, except for the change 
in the X-ray surface brightness which is lowered significantly in the constant energy scaling. 

3.9 Summary of the Analytical Section 

According to our analytic approximations, an explosion of 10^^ ergs into a medium with scale 
height of 22 pc and a density of approximately 6 cm^"^ at the explosion site, when viewed at an 
age of 20,000 years has the following successes as a model of W44: 



1. Its radius is consistent with the angular size and distance. 

2. It will have formed a neutral hydrogen shell on one side (left rear with our orientation), with 
approximately the observed column density. 

3. The expansion velocity of the HI shell will be close to that observed. 

4. With thermal conduction at roughly its unquenched level, the central temperature and density 
are in agreement with those found from the x-ray emission. The dependence of the values on 
the actual magnitude of the conductivity is weak. 

5. The temperature near the edge will be too low to contribute substantially to the x-ray 
emission, so that spectral softening with radius, particularly when combined with the high 
observed column density of intervening absorbers, can provide the appearance of center filled 
x-ray emission. 

6. The van dor Laan mechanism — synchrotron enhancement by compression of ambient cosmic 
ray electrons and magnetic field — can provide the observed radio continuum emission, but 
only with about a factor of 4 higher ambient CR electron population than estimated for the 
solar neighborhood. Additional complicating factors cannot be ruled out. and we discussed 
the possibility of diffusive acceleration boosting the energies of the ambient cosmic ray elec- 
trons in the gas subshock or further acceleration in the converging flows bounding the dense 
shell. 

7. On-going shell formation and the radiative shock on the denser half lead to roughly 2000 
Lq of high contrast filamentary Ha from the dense fast moving shell, consistent with that 
measured after allowing for extinction, plus about 800 Lq of Ha from diffuse fluorescence 
in the surrounding 3 pc or so of stationary ambient gas. The shell radiation is predicted to 
include also the full spectrum characteristic of a fast radiative shock, including amounts of 
[SII] and 63 fiia [Ol] emission completely consistent with those actually observed. 

8. The gamma ray emission expected from the dense shell is 25 to 50% of that seen in the 
vicinity with EGRET. 

4 Physics Assumptions, Justifications, and Implementa- 
tions in the Hydrocodes 

This section surveys the assumptions and methods of the hydrocodes employed in the next two 
sections to examine the spatial configuration of the remnant model's characteristics, its appearances 
in 21 cm, radio synchrotron, and x-ray emission, and its x-ray spectrum. Particular attention is 
paid to justification of the inclusion of thermal conduction in the hot interior of the remnant. 

4.1 Pressure 

The thermal pressure is that of an ideal monoatomic gas with solar composition, assuming full 
ionization of hydrogen and helium, and equal temperatures for electrons and ions. Nonthermal 
pressure is represented via an isotropic term proportional to nP, whose principle effect is to limit 
the compression of the cold dense shell to a realistic value. Although this is referred to as magnetic 
pressure and the density dependence justified by flux freezing and compression of the tangential 
component, the codes do not perform magnetohydrodynamic calculations. Cosmic ray acceleration 
by the shock is assumed to have a negligible effect on the shock or remnant structure. 

4.2 Thermal Radiation 

In the 1-d code, nonequilibrium ionization and the corresponding cooling rates and spectra are 
followed in each parcel with the Raymond and Smith codes (Raymond and Smith 1977, Raymond 
1995), accelerated with modifications using tables of rate coefficients provided by Edgar (1994). 



The abundance set adopted was that of Grevesse and Anders (1988). Radiative transfer between 
parcels was not included. 

Because the 2d code does not follow the ion states with time, it is necessary to estimate 
the radiation rate with a cooling function. As proposed by Edgar (Edgar and Chevaher, 1986), 
a pseudo non-equilibrium radiative cooling table was adopted. It contains the non-equilibrium 
radiative cooling coefficient versus temperature for an optically thin plasma that was shock heated 
to 10^ K (though the choice of upper temperature has little effect) and then isobarically cooled. 

The model x-ray spectra were also generated with the Raymond and Smith code, based on the 
temperature and non-equilibrium ionization state of each parcel previously calculated by the 1-d 
hydrocode. The spectrum was attenuated using the absorption curves of Morrison & McCammon 
(1983), and then convolved with the ROSAT PSPC and Einstein SSS response functions. 

4.3 Dust 

For the most part, dust is assumed to be non-existent within the hot portions of the remnant. As 
shown by Smith et al. (1996), this assumption has little effect on the inferred cooling timescale 
of the gas parcels. The elemental abundances do, however, affect the spectrum. As a result, we 
discuss the results of a 1-d run with dust destruction included explicitly, using the rate model of 
Smith et al. (1996). 

4.4 Thermal Conduction 

The thermal conduction treatment is similar to that of Slavin and Cox (1992), a single temperature 
model with an electron thermal conductivity which is classical in regions of weak temperature 
gradient but bounded by the saturation limit when the gradient is large. As shown in Cui and Cox 
(1992), by the times of interest here this provides a description which is close to that of models with 
different electron and ion temperatures and their corresponding electron and ion conductivities. 

In the present study, however, a single temperature model is not without dangers: examination 
of the Cui and Cox thermal structures shows that there is disparity at the center of the remnant. In 
their case 4, in which the electron and ion temperatures are followed separately with only coulomb 
collisions to equilibrate them but each with its own conductivity, there is a persistent core of hotter 
ions and lower density (about a factor of 3 in the inner third in radius), compared to the single 
temperature case. Normally, this disparity is unimportant; we do not generally imagine these 
hydrocode models to be useful in the central parts of the remnant where the structure is surely 
affected by the details of the explosion and the pre-supernova stellar evolution. In Cui and Cox, for 
example, it is probably a relic of initializing with a Sedov structure. But, when we are purporting 
to model the central x-ray emission, we cannot help being concerned about this point. 

The next concern derives from the considerable uncertainty regarding how much conduction 
actually occurs in the hot interior of an SNR, relative to standard formulae, and how the answer 
depends on time. Although we regard this as a sufficiently difficult problem that observations of 
objects like W44 will probably be called upon to solve it, there are some points worth making. For 
example, it is important to realize that a very small amount of transport, compared to the classical 
heat flux, will make a very large difference to the central structure of a remnant (section |3.3| ) , even 
if that transport is active only when the remnant is rather young. 

The most common arguments against the presence of significant heat transport in remnants 
have to do with the strong suppression of conductivity transverse to magnetic fields. In 2d MHD 
models, a uniform ambient field is distorted outward in the bubble, in most places providing a 
significant tangential component, guiding what conduction there is to the magnetic caps. Arguing 
that there has been considerable turbulence in the center, or that the field was largely irregular 
in the first place, encourages the idea that all conduction paths along the field will be very long, 
suppressing conduction to a large degree. 

Of course, the self-same turbulence that was imagined to have tangled the field would also have 
intermixed the entropy components, accomplishing largely the same purpose as conduction. At the 
very least, it would have smeared a tangled mixture of different temperature components across 
the interior structure. In addition, field line wandering (Jokipii and Parker, 1969) then allows 
transport among seemingly isolated regions. 



In very young remnants, still interacting with their own wind material, there is a lot of radial 
mixing and the dominant field component is radial. This would seem not to suppress conduction 
substantially. On the other hand, for W44, entropy mixing among roughly the innermost 30-100 
M© is needed. 

In point of fact, all conceptions of thermal conduction in this environment are vastly too sim- 
plistic: the collision mean free paths are very long; the stellar ejecta are surely an important 
complication; nonlaminar flow with mixing and turbulence has almost surely occurred; and the 
strict adiabaticity of the Sedov solution is surely the worst possible approximation. Our use of 
classical conductivity may seem naive, but, when safeguarded for saturation, it provides a reason- 
able, physically based estimate for the rate of transport. It could, in fact be an underestimate if 
turbulent diffusion is important in remnant interiors. 

Two papers which have discussed magnetic suppression of conduction and limits on that sup- 
pression include Rosner and Tucker (1989), and Tao (1995). 

In summary, no existing model offers a better idea. This one is important to explore, but 
such exploration must be followed by careful discussion of the discrepancies so we can have some 
indication of what other processes might be important. 
Transport of ions 

A related issue is that in remnant interiors the transport of energy is very likely accompanied 
by the convection or diffusion of elemental abundances and of the ionization structure. This too 
could alter the x-ray spectrum and surface brightness of the center, presumably moving toward 
greater uniformity. 
Implementation of conduction 

In the 2d hydrocode, expressions for the classical conductivity, the saturation flux, and the inter- 
polation between them were taken from Slavin and Cox (1992). The conduction equation is solved 
at the beginning of each base-grid time step on a uniform grid with the resolution corresponding 
to the resolution of the finest grid used in the simulation. An implicit, fractional-steps method of 
Yanenko (1971) is applied. The implementation is similar in Id, except that the originally proposed 
formula for interpolation between classical and saturated fiuxes (-Ftotai = ^ciass/[l + |-Fciass|/^sat]) 
of Cowie and McKee (1977) was used. E xcept as noted, the classical conductivity assumed was 



the full unquenched value as given in § O 



4.5 Synchrotron Emissivity for Use With Hydrocode Results 



Estimation via the methods of Section |3.4.1 of the radio synchrotron emission expected due to 



the van der Laan mechanism requires somewhat unusual accounting, evaluating for each parcel in 
the model its initial volume, compression factor, and gain. The specifics are: the emissivity of the 
ambient medium at 100 MHz was taken as 10"^" W/ (m^ Hz sr); the compression factor was taken 
as X — p/po (where po is the initial density of the parcel); the ambient emissivity was multiplied 
by the initial volume of the parcel, taken to be M/po, dividing by 10® to accommodate the mixed 
units; that initial emission rate was multiplied by the gain at a given frequency, G{x) « (7r/4)a;^ 
for an effective spectral index of 0.5, and then An to eliminate steradians. The final result is the 
parcel's spectral emission (in Watts Hz^^). 

To obtain the total spectral luminosity, the emission is summed over the parcels in the remnant. 
To turn this into a contribution to the observed flux density, it is divided by And? with the distance 
d in meters, to get Watts (m^ Hz)^^ , and divided by 10~^® to convert to Janskys. 



5 Two-Dimensional Model 
5.1 Background 

Purpose 

The 2d modeling provides a fascinating look at the structure of the remnant as well as quan- 
titative representations of the HI distribution, across the remnant and in radial velocity, and of 
the radio continuum emission distribution anticipated from the van der Laan mechanism. It also 
allows us to examine the x-ray emission from the distribution of density and temperature in the 



interior to see whether lack of sphericity might invahdate our more detailed Id model of the x-ray 
surface brightness distribution and spectrum. 

Code Description 

The hydrodynamic and thermal conduction equations are solved with the help of the AMRA 
code (Plewa & Miiller 1998), which combines the AMR (adaptive mesh refinement) algorithm 
of Berger and Colella (1989) and the PPM method of Colella & Woodward (1984). It has also 
been recently used in a numerical study of the evolution of ejecta fragments by Cid-Fernandes 
et al. (1996). The specifics of our implementation of thermal conduction, nonthermal pressure, 
and thermal radiation are discussed in section ^. (The effective magnetic field strength assumed 
was quite low, corresponding to 2.1 /iG at 20 cm~'^ and proportional to density, so the nonthermal 
pressure was felt only in the dense shell and not very strongly even there.) 

The present application employs spherical coordinates (r, 0, </>) and rotational symmetry with 
respect to the polar {9 — 0, 180°) axis. The initial grid extends over radii 0.0 < r < rmax and polar 
angles 0° < 6* < 180°, containing 128 and 15 uniformly distributed points in r and 9, respectively. 
The adaptive mesh refinement takes place in two stages, as needed, with total factors of up to 20 
and 6 improvement in r and 9 resolution. Thus, the resolution of the finest grid is equivalent to 
that of a uniform grid of 2560x90 points extending over the whole computational domain, while 
requiring about a factor of 15 less CPU time. This resolution was chosen after Id experiments 
showed that with 2560 zones in radius the fully formed shell will be resolved into not less than 10 
cells. At lower resolution we would not have been be able to describe the velocity distribution of 
the cold material. Even with the adaptive mesh, the CPU time used for the simulation presented 
below amounts to ~ 10 days of the CRAY J90. 
Initial conditions 

The code is initialized at a remnant age of 300 years, with the explosion energy in expansion 
of the ejecta. The mass of the ejecta is 5 M© ; the density and velocity profiles are those used by 
Cid- Fernandes et al. (1996). The outer expansion velocity is about 5000 km s^^ and therefore 
the radius of the ejecta at initialization is about 1.5 pc. 

The values assumed for the explosion energy (10^^ ergs), ambient density at the explosion site 
(6.2 cm"'^), and scale height of the ambient density distribution (22 pc) are those found from our 
analytical study as recorded in Table 

5.2 Results 

Each cell, or parcel, consists of an annulus in 0, with volume 27rr^ sin 9drS9. At a specific evolution- 
ary time the output for each annulus consists of the average r and 9, along with the corresponding 
density, temperature, mass, and the two components of velocity. From these, the thermal and mag- 
netic pressures can be inferred, along with the x-ray emission (assuming coUisional equilibrium) 
and radio continuum emissivities. 

5.2.1 Structure 

Figures ^, ^, and ^ show the density, temperature, and pressure structures of the simulated remnant 
at an age of 20,000 years. In the transparent side view of Fig the densities in all cells are shown, 
versus their locations along the major axis. The locus of values, between 3 and 10 cm~"^, making up 
the straight diagonal line extending beyond the remnant shows the exponential density distribution 
of the ambient medium. The diameter (major axis) is seen to be just over 22 pc, essentially one 
scale height in the ambient medium. 

There are deviations from spherical symmetry in the low density interior, but they are quite 
modest. The central density is about 15% of ambient at the equator, its minimum being about 1.0 
cm~^, as predicted by the analytical modeling. In addition the distribution is very flat, exceeding 
the minimum by less than a factor of two over a large region in the interior. 

The cells with very high density are those in the cold dense shell. As we discuss below, the 
thermal pressure was too large at high density in this run, and the magnetic pressure was too low. 
With these conditions, the densest ones should have been compressed by about a factor of 120 over 
ambient, and are, as shown. Complete shell formation, shown by the locus of highest density cells. 



has reached from the dense end to the equator, just as predicted, beyond which the maximum 
compression fades sharply. 

The many cehs with compressions somewhat less than maximum show the well resolved tran- 
sition of material which has begun cooling but has not yet reached its limiting density. At the 
tenuous end, shock compression by only the adiabatic factor of 4 is apparent, while not far away 
where radiation is just beginning to affect the structure, that factor of 4 is followed by further 
compression. The final distribution of grid resolution is also apparent in the figure, where the 
regions of highest density have the greatest concentration of cells in r and 9. 

The temperature structure in Figure || is shown as a 3d surface plot with the densest end 
on the near left side. It shows a central temperature of about 6 x 10^ K, consistent with our 
expectation from the conduction analysis (Section 3.3), with a slightly lower shoulder just left 
of center where the density was higher. The trough in temperature behind the shock on the 
tenuous end is a surprise, indicating that the shock there is already becoming radiative where our 
analysis had predicted a cooling time of 25,500 years (Table ||). This is consistent, however, with 
indications (below) of a relatively low pressure on the tenuous end. As discussed in Maciejewski 
and Cox (1998), a pressure difference between the two ends is apparently also responsible for the 
reduced displacement observed for the remnant center. Apart from the aforementioned shoulder 
and the radiative shock girding the tenuous end, the temperature structure is nearly symmetric 
and parabolic. Both the central value and the general distribution are quite similar to Id results 
presented in Section ^ below. 

The pressure in the interior of the remnant (Figure ^ is very flat in the center with a value of 
1.6 X 1G~^ dyn cm~^, but drops abruptly by 5% about halfway to the edge. On the tenuous end, 
this drop leads to a broad shallow valley at 85% of the central pressure. On the dense end the 
pressure continues to drop with increasing slope, to a deep trough. 

Very close to the edge on the dense end, it is apparent that the low pressure is due to rapid 
radiative cooling. The structure further in, which can also be discerned in the density details 
of Figure ^ resembles a substantial rarefaction moving inward or a conipressional wave moving 
outward. Such waves are fairly commonly found in hydrodynamic models although none is obvious 
in the particular Id model presented in the next section. When present, they have very long 
reverberation periods and thus appear for a long time as nearly stable features. We explored several 
possible sources for this structure in our 2d run and found that the high flat central pressure, its 
weak but abrupt edge, and its asphericity are the residual of the reflected reverse shock. (An 
experimental Id run with the same initialization as the 2d run produced a very similar structure.) 
Because a reverse shock is required to raise the pressure of the center, and it will surely be reflected 
to some degree, there should be a substantial time during which this sort of structure is present 
within a remnant. But detailed conclusions based on its presence in a particular remnant model 
should be regarded with some suspicion. When it is found, where, and with what strength, depends 
on the initialization. Perhaps a fruitful approach would be to search the x-ray emission pattern of 
a remnant for indications of such a structure and, if found, to use the structure to evaluate models 
of the remnant's early state. 

Stating our pressure results more optimistically, the interior pressures in our standard Id and 
2d runs differ by less than 30% anywhere that the temperature exceeds 2 x 10^ K. The differences 
will have negligible impact on the interpretation of the x- ray data. 

The pressure structure at the very edge of the remnant bears an imprint of both the cooling 
history and the way we treat cooling in 2d. (Note that the pressure spikes shown are an artifact 
of the way in which the figure was prepared, not a characteristic of the model. They are due to 
incomplete sampling of the densest region in converting the model data from spherical to cartesian 
coordinates for the plot. This makes the shell conveniently transparent, but inaccurately suggests 
that it was poorly resolved in the model.) At the tenuous end, the pressure maximum is in the 
immediate post shock gas, prior to the onset of cooling. Further inward where cooling has begun, 
the pressure drops precipitously. The structure at this end is entirely resolved. At the high density 
end, the radiative shock is unresolved and the high pressure region is in the already cold dense 
shell, itself trailed by a deep pressure trough of cold but not yet compressed material. Comparing 
the maximum pressures at the two ends, we found 3.4 x 10^^ and 2.4 x 10~^ dyn cm~^, for the 
dense and tenuous ends, respectively. Moving away from either end the edge pressure drops rather 
rapidly. On the tenuous end this appears to be due to the onset of cooling, which coupled with the 



inverse temperature dependence of the cooling coefficient, temporarily stalls the shock, possibly 
leading to a period of instabihty and secondary shock formation (Falle, 1975, 1981). This is most 
exaggerated in an extended region of the remnant's surface over which the post shock pressure 
seems to be very low, just toward the tenuous end from the equator. In this region the radiative 
shock appears not to be resolved (the initial compression exceeds 4), but cooling is so rapid that 
the post shock pressure is small and the shock velocity temporarily low. On the dense half, the 
pressure in the cold shell drops a factor of two between the end and the equator. On the whole, 
these pecularities in the pressure structure provide the biggest deviations from the expectations 
based on analytic and Id models. 

In summary, the magnitudes of central density, temperature and pressure are all as expected. 
Similarly, shell formation over half the surface has occurred as planned. There were some new 
but qualitatively understandable features brought to light, particularly the strong variation of 
edge pressure over the surface at this epoch, and the seemingly premature onset of cooling on the 
tenuous end. We are left not knowing whether the detailed interior pressure distribution of the 2d 



case is a better or worse approximation than the Id structure presented in Section 6.2.1, but are 
relieved that the differences are modest. 

5.2.2 HI Shell 

In the particular orientation of our 2d model, the symmetry axis is tilted 40 degrees from the line of 
sight with the dense end on the far side, and toward the northeast, giving us a glancing view into the 
open end of the half shell. With this projection, most of its dense material (and associated strong 
magnetic field) is in the northeast, where the average radio continuum brightness is highest, and 
the dense gas is moving predominately away, so that the ranges of expansion velocities observable 
via 21 centimeter emission are occupied only on the receding side, as observed. 

In preparing the figures of this section, the mass of each r, 9 annulus of the calculation is 
represented by a proportional number of dots, randomly distributed within the annular volume. 
The spatial coordinates and velocity components of each dot are then tracked through the rotations 
to the final projection. (In later sections on the radio continuum and x-ray appearances, the process 
is similar but the number of dots per annulus is chosen in proportion to emission rather than mass.) 
The dot density in the shading thus represents the projected surface density (or surface brightness 
in the emission pictures). Unless otherwise noted, the assumed distance of W44 was that required 
to give an angular diameter of 35', about 2200 pc, closer even than our revised distance estimate 
as quoted in Table |l|. Methods for scaling the results were presented in Section 3.8 and their 
application will be discussed further in the conclusions. 

Figures |^ and |^ present our results for the distributions of mass and velocity in the dense shell, 
in formats for direct comparison with the observations of Koo and Heiles (1995). They include only 
cold material more dense than 100 cm~^, assuring a sufficiently short (^ 1000 year) recombination 
time to achieve neutrality (the ionization state not actually having been tracked in the model). 

Figure ^ consists of 4 panels, two of which are projected column density, the others being 
position velocity diagrams. The upper right panel shows the total column density map of the 
dense half shell. As expected, it is concentrated toward the northeast, sharply bounded and edge 
brightened there, but fades out gradually to the southwest. 

This view, however, is not available to observers because of the large amount of low velocity 
hydrogen toward W44. The corresponding view presented by Koo and Heiles (1995) is limited to 
the velocity range i^lsr = 135 to 170 km s^^ and displayed as their Figure la. Their figure and the 
corresponding model column density map for that velocity range (including the systemic velocity of 
+43 km s~^) are shown in the lower left (oblique view) panel of Figure ^ The model view changes 
considerably when limited to this velocity range. It is much more diffuse, lacking the sharp bright 
northeastern edge, smaller, more symmetric, and has a hole in the center from omission of highest 
velocities. With these changes, it bears a very strong resemblance to the HI observations. 

The model's restricted velocity view has a total diameter of about 25' (a mass weighted diameter 
of about 18') versus the shell's total diameter of about 35'. For comparison, Koo and Heiles quote 
the mass weighted radius at wlsr = 140 km s~^ as 7.5'. The corresponding mass weighted 15' 
diameter is smaller than our model. As will become apparent below, the model is in good agreement 
with the data envelope in position velocity diagrams, but the HI data are more filled in velocity 



space, leading to a lower mass weighted size. 

Taking model data from within the apertures outlined in the upper right panel of Figure |, the 
upper left and lower right panels show the corresponding position velocity diagrams. In contrast 
to the total surface density plots, the half shell of the model is immediately apparent. The model 
remnant would be observable in 21 cm almost exclusively on the rear (receding) surface, as is W44. 
(If our model were absolutely correct in its degree of shell formation and orientation, however, we 
would predict the observability of a small quantity of rapidly approaching HI in the northeast.) 

In Figure ^ portions of the position velocity diagrams with W44's systemic velocity of +43 
km/s added, are replotted with the same restricted velocity range as in Figs. Ic and Id of Koo 
& Heiles (1995). In addition, the HI data were kindly supplied to us by Koo, and are shown for 
comparison. 

The correspondence between the observations and the model is again very good: both show 
the highest velocity of the bulk of the material to be about 180 km/s, but with some moving as 
rapidly as 200 km/s. The highest velocities arc found toward the center, with lower values near 
the edge and a more hollow center, indicative of an expanding shell. The radial scales correspond 
very well on the whole, with the model distribution outlining the data and having a diameter of 
about 25' at 140 km s~^, about 2/3rds of the actual remnant diameter. 

As noted above, the observations seem to have a wider range of velocity at a given radius than 
does the model, with a lumpier distribution of material. The model is highly idealized, neglecting 
small scale irregularities in the ambient medium and any turbulence that may be present. As 
we discuss and demonstrate below, the radio continuum maps show strong filamentation that is 
associated with surface irregularities which, in addition to their effect on the surface brightness 
profile, will tend to broaden and disturb the velocity structure. For example, a bright radio 
continuum filament indicates a tangential view of that portion of the remnant surface. When such 
a filament is seen projected well away from the remnant edge, the line of sight velocity there will be 
very small, because the surface motion is turned crossways. Any remnant whose radio continuum 
map shows an intricately textured surface like that of W44 will have a much broader and more 
disturbed radial velocity distribution in 21 cm than any smooth shell model. 

Considering that the model made no attempt to include any irregularities, it is somewhat 
surprising that its velocity width at the remnant center is as broad as it is. This appears to be 
due to the ongoing rapid development of the shell. There is material which has cooled and been 
compressed, but has not yet merged with the densest gas in velocity space. This is indicative of the 
strong velocity convergence in the neighborhood of the densest regions, as discussed in Section 3.4.2| 



with regard to suppression of cosmic ray diffusion from the shell. There may also be a contribution 
from the instability of radiative shock waves at high velocity (Fallc, 1975, 1981). 

Koo and Heiles (1995) reported that the mass associated with the shell and having wlsr > 130 
km s~^ is 73 ± 3 M0. The mass of the simulated shell having velocities within this range is 149 
Mq at an age of 20,000 years. The dependence of the model results on remnant age are shown in 
Table | both for the mass within that velocity range and for the total mass at density above 100 
cm~^. (The results depend only slightly on the density criterion for inclusion in the neutral shell.) 
Apparently the model shell's cold mass is about a factor of two larger than the observed HI mass, 
a conclusion that does not depend strongly on the age at which it is evaluated. 

Several neglected factors are probably responsible for this differential. For example, the shell 



is partly ionized. In Section 3.5 we found that (because of photoionization) each atom joining the 
shell had to recombinc about 3 times. We expect the ionized gas in the shell to be roughly equal 
to that collected during the past 3000 years (3 recombination times at 100 cm~'^), lowering our 
estimate of the neutral mass with wlsr ^ 130 km s^^ by up to 50 Mq . Similarly, the presence of 
the OH masers implies that the shell is at least partly molecular. If 200 of the 450 M© in the shell 
were H2, there would be no discrepancy. Finally, the fact that the observed remnant has surface 
structure that broadens the velocity distribution in its position velocity diagrams makes it difficult 
to compare the amounts of mass in an observed velocity range. The distortion could easily have 
skewed the velocity distribution downward near the cutoff velocity, reducing the measured shell 



Table 3: Time Dependence of 2D Model Characteristics 



Age (years) 


16,000 


17,000 


19,000 


20,000 


21,000 


Cooled Mass (Mq ) 


230 


310 


475 


565 


655 


Mass in Shell with 












n > 100/cc (Mq ) 


135 


210 


370 


450 


535 


Shell Mass with 












V(LbRj > 130 kni/s (M© ) 


60 


88 


134 


149 


157 


Radio Continuum 


104 


179 


346 


423 


503 


100 MHz (Janskys) 












ROSAT Hard Band (counts/s) 












NH = 1.0E22 


8.5 


7.3 


5.5 


4.8 


4.2 


NH = 1.5E22 


4.8 


4.1 


3.1 


2.7 


2.4 


NH = 2.0E22 


2.9 


2.5 


1.9 


1.6 


1.4 


ROSAT Soft Band (counts/s) 












NH = 1.0E22 


29.0 


24.9 


18.9 


16.6 


14.7 


NH = 1.5E22 


7.5 


6.5 


4.9 


4.3 


3.8 


NH = 2.0E22 


2.5 


2.1 


1.6 


1.4 


1.2 



5.2.3 Radio Continuum 

This section discusses the radio continuum emission expected due to the van der Laan mechanism, 



following the prescription in Section 4^. Before doing so, however, we had to overcome one 
difficulty. Our 2d model assumed an unrealistically low value for the magnetic pressure. In addition, 
cooling was suppressed below 10"* K, yielding an excessive thermal pressure. Together these led to 
results for density in the shell that were unreliable (though the two effects roughly cancelled and 
the density structure turned out to be reasonable. For this section, we corrected the compression 
factor in each annulus of the calculation, under the assumptions that: the hydrocode provides 
a correct value for the annulus's total pressure, the initial density is the ambient density just 
outside the shell at the same z (that is, that material in the dense shell has not moved far in 
z), the ambient tangential field outside the remnant is 3/iG independent of external density, and 
that thermal and cosmic ray pressures are negligible within the shell (yielding a reasonable but 
optimistic assessment). Setting the magnetic pressure equal to the total computed pressure in the 
annulus determines the corrected compression factor for calculating the gain. Following the rest of 
the prescription provided the flux densities shown in Table |^. As found in the analytic section, this 
slightly optimistic estimate of the van der Laan flux density is consistent with the observation of 
500 Janskys at 100 MHz. But once again there is the hidden assumption that the background can 
be produced by CR electrons in the mean ambient field, an assumption that fails by approximately 
a factor of 4 when applied to estimates in the solar neighborhood. 

The 2d model also makes possible a projected image of the synchrotron emission, prepared 
in much the same way as was done for the HI structure. Under the assumption that the true 
radio continuum emissivity is proportional to our van der Laan modeling, and neglecting the 
directionality and polarization of the emission relative to the magnetic field, our map for the radio 
continuum is shown for our standard projection in the upper right panel of Figure |^. A similar 
image in the lower left panel shows the view that would be had from a direction perpendicular to 
the major axis; the half shell of the model is obvious. The figure also shows three horizontal and 
three vertical slices through the emission. The central slices quantify the strong degree of edge 
brightening apparent on the northeast. 

The observed radio continuum distribution, illustrated beautifully in Smith et al. (1993) and in 
Giacani et al. (1997), is qualitatively similar but differs in two significant characteristics. The ob- 
servations in fact were the stimulus for our projected half shell model, giving the distinct impression 
of a view into an open hemisphere, brighter in the northeast than southwest because of the double 
contribution from the two sides in projection. Our model, however, has a blurred edge extending 



further around the sheU, and thus projects into a more rounded configuration with less of a step 
in overall brightness. The second major difference is that the observations show quite pronounced 
filamentation, some of which appears to be associated with the sharpness of the equatorial edge. 
Hester (1987) has shown in the case of the optical emission of radiative shocks that a high degree 
of filamentation can derive from spatial irregularities in very thin emitting structure. In order to 
test this as a possible source for the filamentation of the radio continuum in W44, we prepared 
an image of the 2d model in which the individual emission elements were radially displaced in a 
sinusoidal pattern resembling the corrugated surface of an acorn squash. The results are shown 
in Figure ^ and very clearly show that a modest degree of surface crinkling can lead to striking 
filamentation. The horizontal and vertical slices show the filamentation contrast. 

In modeling the filamentation in this way, the areas of especially bright radio continuum are 
assumed to derive from tangencies along the surface from regions which are extremely thin. As 
such, they should also be filaments in HI, optical, and all other emissions associated with the dense 
shell and radiative shock, as is observed. 

Because the ambient cosmic ray electrons have a flatter energy distribution at lower energies, 
the radio continuum spectrum of W44's filaments should be flatter than that of the remnant's more 
diffuse component. This comes about because, qualitatively, there are two layers of radio emission, 
a very thin bright one from fully compressed material in the dense shell, and a somewhat thicker 
layer from cooled material in the process of compression. The thinner layer dominates the filaments 
because it is most edge brightened by tangency, while the two layers contribute more equally to the 
diffuse component. On average, then, the filaments have suffered greater compression and their 
emission derives from electrons which were initially at lower energy, where the spectrum is flatter. 
Potential applications of this phenomena include the possibility of inferring the very low energy 
interstellar cosmic ray electron energy distribution from filament spectra of very high pressure 
(small, high density) supernova remnants. 

5.2.4 X-Rays 

The x-ray surface brightness distribution and count rate for the 2d model were calculated from 
equilibrium Raymond and Smith spectra, folded through ROSAT PSPC response functions for 
two bands (soft R4+R5+R6, 0.4 - 1.56 keV, and hard R7, 1.05- 2.04 keV), and attenuated by 
photoelectric absorption by the foreground material. The total count rates for both bands, at 
three absorption column densities, for five ages are shown in Table ^. The observed band ratio of 
about 1 for these two bands (Rho et ai, 1994) points to an absorption column between 1.8 and 1.9 
xlO^^ cm~^. At an age of 20,000 years, the corresponding total predicted ROSAT PSPC count 
rate is then also close to the observed 4.2 s~^ found by Harms et al. (1997). 

The distribution of X-ray surface brightness for this case is shown in Figure ^. It is brightest 
in the center, approximately gaussian in appearance, and highly symmetric about the center. The 
full width at half intensity is about 19', just over half the full diameter. There is nothing about 
the structure that appears to be particular to the two dimensionality of the calculation or lead 
one to suspect there was a significant density gradient in the ambient medium. The center of the 
emission is shifted slightly away from the explosions site, but is well centered within the remnant. 
Of course, when comparing with the observations, one has to bear in mind the strong sensitivity 
of the count rates to the absorption column density. It is quite possible that variability in the 
absorption screen leads to artificial structure in the observed x-ray surface brightness. Scattering 
of the X-rays by intervening dust may also be a factor. 

The appearance the model remnant would have without absorption is shown in Figure |l^. 
From this result, it is clear that at this particular age, the outer parts are everywhere too cool 
to contribute to the emission and the absence of the outer shell emission is not primarily due 
to absorption. The fact that the pressure drops somewhat away from the remnant center also 
contributes to this behavior. 

^From these results we are confident that ID modeling of the noncquilibrium ionization struc- 
ture and x-ray spectrum presented in the next section reasonably corresponds to a hypothetical 2d 
model which followed the ionization. 



Table 4: Parameters of Id Runs 



Uo = 6.2 atoms cm 


£■0 = 1.0 xlO^^ ergs, 


t = 20,000 years 


Model 


Conductivity 


Dust Destruction 


Abundances 


A 


full 


utter 


solar 


B 


zero 


utter 


solar 


C 


full 


utter 


enhanced 


D 


zero 


utter 


enhanced 


E 


l/lOth 


utter 


solar 


F 


full 


sputtering 


solar 



6 One-Dimensional Model 
6.1 Background 

Purpose 

A variety of Id hydrocode models were made during the course of this investigation. They were 
used to verify the approximate cooling time formula, to explore the lag time between cooling and 
shell formation described in Section |3.2| , to explore the time development of shell mass and radio 
continuum emission, to measure the resolution needed in the 2d model, to examine the dependence 
of the central conditions (density and temperature) on the assumed degree of suppression of the 
conductivity, to explore the degree of dust destruction expected within the hot interior of the 
remnant, to simulate inclusion of additional heavy element abundances in the ejecta, to explore the 
dependence of the interior pressure structure on the SNR initialization characteristics as per Section 



5.2.1, and to calculate non-equilibrium x-ray spectra for comparison with the surface brightness 
distribution and with the available spectral information. 
Code Description 

Our code architecture employs a Lagrangian mesh, based on the scheme discussed in Richtmyer 
and Morton (1967) and described in detail in Smith and Cox (1998). The method is similar to that 
of Cui and Cox (1992) and Slavin and Cox (1992). The ionization and emission characteristics are 
discussed in Section 
Initial conditions 

A simulated supernova explosion is initialized as a 2 pc radius sphere containing 2 solar masses 
of ejecta with a velocity linearly proportional to the distance from the center. 

The values assumed for the explosion energy (10^^ ergs), and ambient density at the explosion 
site (6.2 cm^'^) are those found in our analytical survey and recorded in Table ||. All results 
reported are for an age of 20,000 years. With these parameters and age, the model SNR is found in 
an adolescent stage of shell formation. It approximates edge conditions only near the equator of the 
2d model, and is therefore not directly useful in studying the radio continuum or shell development. 
We have, however, examined these characteristics at later times and in models that were higher in 
density to verify behaviors of the 2d results. 

Results of six models appear below, differing in assumptions about the level of thermal con- 
ductivity (full, reduced, zero), dust destruction (utter, thermal spallation only) and abundances in 
the ejecta (solar, 5 solar). Their different characteristics are summarized in Table ^. 

6.2 Results 

At a specific evolutionary time the output for each parcel consists of the inner and outer radii and 
velocities of each parcel, plus averages of density, temperature, thermal and nonthermal pressures, 
and nonequilibrium ionization. In models following gradual dust destruction, or enrichment in 
ejecta, the gas phase elemental abundances are also recorded for each parcel. From these, the 
x-ray spectral emissivity of the hot interior, shell mass and its velocity distribution, and radio 



continuum emission can be calculated as described in Sections 4.2 and 4.5 



6.2.1 Structure 



Figures |T^, and |l^ show the density, temperature, and thermal pressure structures of Models 
A and B at an age of 20,000 years. The differences, due entirely to thermal conduction, are 
immediately apparent.^ In agreement with the 2d run, the radius is just over 11 pc and the 
interior electron density in Model A is about 0.92 cm~^ (nuclear density about 0.84 cm~'^), and 
very flat. This central density is 13% of the ambient density, just as predicted in Section |3^. For 
Model E (not shown), with its thermal conductivity reduced by a factor of 10, the central electron 
density is still fairly flat but has a value of only 0.51 cm~^. The central density reduction by a 
factor of 1.8 with a conductivity reduction of 10 is also exactly that expected from the analytical 
work in Section 3.3. 

The temperature structure in the interior of Model A also reconfirms the 2d model, dropping 
only gradually with radius from a central value of 6x10^ K, in the range consistent with x-ray 
observations. In the outer parts, thermal conduction has delayed the progress of shell formation 
somewhat and provided a slightly larger remnant. Models A and B both have similar temperature 
and density ranges in the outer parts, so that, as in Harrus et al. (1997), it is possible to find 
an x-ray model with a similar total luminosity, but it will be flatter or even shell like in surface 
brightness, and may exhibit a stronger radial temperature profile. 

The pressure in the remnant interior is very flat in either case, as expecte d, wi th a value of 
about 1.5xl0~^ dyn cm~^ for Model A. As previously mentioned in Section ^.2.1 the pressure 
distribution of Model A is smoother than that of the 2d model and the difference was traced to 
the different initializations. 



6.2.2 X-ray Surface Brightness Distribution 

Figures and b compare the predicted radial flux distributions from Models A and B (after 
absorption) with observations in the ROSAT PSPC soft (R4+R5+R6, 0.4 - 1.56 keV) and hard 
(R7, 1.05- 2.04 keV) bands. The model spectra were attenuated by a constant Nh = 1.89 x 10^^ 
cm~^ (found by fitting the Model A spectra to observations, below), although Rho et al. (1994) 
found best fit values for coUisional equilibrium isothermal models varying between 1.57 and 2.10 
xlO^^ cm~^ across the remnant. The model remnant was taken to be at a distance of 2.2 kpc 
for these plots so that the angular diameter of the shock would be 35'. This has no effect on the 
magnitude of the surface brightness but determines the angular scale and PSPC count rate. 

Obtaining a meaningful radial distribution from the asymmetric ROSAT image involved several 
steps. The data were cleaned and calibrated using the method described in Snowden et al. (1994), 
and following the examples given in Snowden (1994). Emission from a bright G star on the left 
side of image was then removed and, using Jones et aZ.'s (1993) choice for the nominal center 
{a, 5 = 18''53'"30*, 1°15 (1950)) the data were averaged over a 90° angle in each of the four 
cardinal directions (north, south, east, west). The persistence in Figures [l4| a and b of the surface 
brightness to greater radii in the north is real: the X-ray distribution is longer north-south than 
east-west and the point chosen for the nominal center is south of its centroid. 

Overall, the count rate and surface brightness distributions of Model A and our 2d model are in 
excellent agreement. Compared to observations, they account for several key features: the surface 
brightness is highest at the center, the limb is not brightened, and the total count rate of Model A is 
just 12% higher than observed (see Table |^). Model B, however, without conduction, has a factor of 
2 too little emission overall and is edge brightened. (A model similar to B, by Harrus et al. (1997) 
was somewhat more successful.) The only apparent failings of the conductive models are that their 
emission distribution is dome-like while the observations resemble a peaked volcanic conej^ and 
that the 2d model is too symmetric, not providing the north — south elongation apparent in the 
image. The next few paragraphs record various experiments which explored model complications 
that might sharpen the central peak of the distribution. 
Emission from Thermalized Ejecta? 

^ The glitches at a radius of 6 pc in the Model B results are at the contact discontinuity for that 
run and arise from the discontinuity in velocity at initialization. For Model A, the 2 Mq of ejecta is 
all interior to 2.5 pc and the discontinuity has been erased by thermal conduction. 

^ Note that sharpness of the central brightness peak varies among this class of supernova remnants, 
with W44 and W28 having sharper peaks than 3C391 or W63 (Rho, 1997). 



At 1 crn"'^, a radius of 4 pc contains 8 M0 . Thus, the mass in the central few parsecs of W44 
should clearly be heavily contaminated by the stellar ejecta and enriched in metallicity. Instabilities 
and turbulent mixing potentially offer the possibility of diluted enrichment through much of the 
roughly 100 Mq in the hot interior. With an enrichment that decreased with radius, the observed 
radial profile could be matched. What is surprising, however, is how little enrichment could be 
tolerated in conductive models before the prediction exceeds the observed brightness. For example, 
the factor of only 5 enhancement over solar abundances in the 2 Mq of ejecta in Model C increases 
the line emissivity by a factor of 5 in the inner 2.5 pc ^4', and the central surface brightness by 
about a factor of 2. More than that would be too much unless compensated by modeling a less 
dense remnant at greater distance or by a reduction of thermal conductivity. Interestingly, the 
nonconductive models do not suffer from this limitation; they are so hot in the interior that their 
(already weak) central emission is mostly due to bremsstrahlung and does not increase substantially. 
Our tentative conclusion is that the elemental abundances of the ejecta must be accommodated 
much more realistically before the very centermost region can be addressed. A major part of 
getting that right will be to manage to hide most of the metallicity, in clumps that were never fully 
heated, in overdense regions which have already managed to cool, or in dust grains. 
Cloudlet Evaporation 

Setting aside concerns over the ejecta, we imagined a hybrid model, including thermal conduc- 
tion in the interior (represented by our hydrocode results) and cloudlet evaporation as parame- 
terized by White and Long (1991). In particular, we looked for possibilities in which evaporation 
contributes to the brightness only in the central few parsecs. We expected that somewhat lower 
mass loss from evaporating clouds would be required than estimated by Rho et al. (1994) or Jones 
et al. (1993). This possibility is not clearly successful, however. Figure 4 of White and Long (1991) 
shows that the effect^ of an evaporating population of clouds is to increase the column emission 
measure for projected radii less than about 0.65 to 0.85 of the shock radius (depending upon the 
chosen evaporation parameters). Only in the extreme case of an infinite evaporation time scale 
and very large ratio of preshock cloud mass to preshock intercloud mass is the contribution in the 
center significantly greater than that at medium impact parameters, and even this extreme case is 
not significantly more centrally peaked than our models, owing to their flat central distributions 
of density and pressure. Their model as it stands does not seem to contribute usefully to peaking 
of the distribution. 

Restriction of Dust Destruction to Sputtering 

Restricting dust destruction to thermal sputtering in the hot gas was considered in Model F, 
using the method described in Smith et al. (1996) to calculate the rate. We found that dust 
destruction was greater in the center of the remnant, providing a relative enhancement of the 
abundances of the hot gas, but that the effect was not strong enough to match the peaking of 
emission in the center of the remnant. Additionally, since the plasma with depleted abundances is 
less efficient at X-ray emission, the total flux was reduced to about 50% of that of Model A with the 
same absorption column density. A more careful study including dust destruction processes in the 
shock and revision of absorption based on the new model spectrum might be more successful overall. 
In any case it is somewhat dangerous to assume that the elemental abundances everywhere in the 
hot interior are undepleted. Thermal sputtering alone is insufficient to destroy the dust completely. 
How Can a Remnant be Shell-Like in Radio and Centrally Bright in X-rays? 

Like W44, our Model A does not appear limb brightened. And, as with our 2d model, this 
absence of limb brightening is not dependent on preferential absorption of softer x-rays emitted by 
the cooler material further out, although that effect certainly shapes the profile somewhat. The 
density rises toward the outside, and there is considerable emission from material there, but not at 
x-ray wavelengths. The reason that our model of W44 is centrally brightened in x-ray emission is 
that the explosion occurred in a sufficiently high density that the X-rays emitted by the interior are 
still bright when the outer part of the remnant becomes too cool to detect. Using our rule of thumb 
that, at the shell formation time, the central density in the presence of thermal conduction should 
be about 0.13 of the ambient density, explosion site densities in excess of a few cm"'^ should lead to 
radiative remnants (shell like in radio continuum) with observable x-ray emission from the interior. 

The environment for their work was a model which included thermal conduction only in the the 
evaporation of clouds, while our picture assumes thermal conduction everywhere. Perhaps the two 
views do not mix as easily as we have assumed in this discussion. 



Without thermal conduction, x-ray emission is confined to an interior thick sheU-hke region, as in 
our Model B, with temperature and density similar to Model A. Harrus et al. (1997) showed that 
such a model would not necessarily show its shell like character in projection, although our Model 
B does. Although Model A works better, it still fails to address the central peak without further 
complications. 

6.2.3 X-ray Spectrum 

X-ray Spectrum 

The observed spectra from the ROSAT PSPC and Einstein SSS are compared with Model A in 
Figure and b. Because the SSS had a 6' diameter field of view, the model SSS spectrum was 
calculated using a cylindrical slice through the center of the remnant while the PSPC spectrum 
used the entire remnant. Because our models are somewhat underluminous in the center and the 
Einstein SSS and ROSAT PSPC have different fields of view, the normalization for each instrument 
fit was allowed to vary independently. The total fit to both datasets has a reduced of 2.0. Most 
of the comes from the Einstein SSS data because the model has too little emission in the Si 
XIII complex around 1.85 keV. The best fit normalizations turned out very well. The model was 
underluminous compared to the Einstein SSS observations, by 36%. For the PSPC, however. Model 
A predicts 25% more emission than is seen. Since the SSS had a smaller field of view than the 
PSPC, this result suggests the central region is brighter than this model predicts. The best fit Nh 
was 1.89 X 10^^ cm~^, a value we've chosen to use for all other comparisons, and which concurs 
with the 2d model conclusions. 

The spectrum for Model C with its enhanced abundance ejecta is shown in Figure |l^. The 
column density fit to the PSPC and SSS data is barely changed (Nr = 1.88 x 10^^ cm~^), but the 
PSPC response to the model now implies the model is 35% too bright. The model also overestimates 
the Einstein SSS emission, predicting 23% more emission than observed. The overall spectral fit to 
the data is slightly worse; the reduced = 2.5, and again the worst part is the fit to the Einstein 
SSS at 1.85 keV. 
Ionization Equilibration 

Rather than a thorough exploration of coUisional ionization equilibration (CIE) of the diffuse 
interior in the presence of substantial conduction, we content ourselves with a comparison of the 
x-ray spectra of Model A, with and without the assumption of CIE, to learn what instrument 
characteristics are required to distinguish between them. Although we have the complete nonequi- 
librium ionization (NEI) balance saved for each parcel, we can calculate the corresponding CIE 
model as well. We did two different model fits. In the first case, we allowed both the overall 
normalization and the absorbing column density Njj to vary. In this case, the best fit value (with 
reduced = 1-89) had Nr = 2.03 x 10^^ cm-^. The PSPC model was 10% underluminous, while 
the Einstein SSS model was underluminous by a factor of 2. Figure ^ compares this equilibrium 
model with the nonequilibrium spectra as seen by the ROSAT PSPC and the Einstein SSS. For 
the second case we held the column density fixed at Nh = 1-89 x 10^^ cm^^ and fit the data. In 
this case the fit was extremely poor (reduced = 4.7). The PSPC model, unsurprisingly, is now 
overluminous by 14% while the Einstein SSS model emission is still a factor of 1.82 underluminous. 

The Einstein SSS spectrum shows three strong peaks, which were identified by Jones et al. 
(1993) to be from Mg XI (1.34 keV), Si xiii (1.85 keV), and S xv (2.45 keV). These ions are all 
part of the helium isosequence, and each corresponds to the strong helium triplet of forbidden, 
intercombination, and resonance lines. (See Gabriel et al. (1991) for a fuller discussion of these 
lines in the context of O vii.) The largest difference between the CIE and NEI models is that in 
equilibrium, the dominant stage is helium-like. In the NEI model, the gas has been over-ionized by 
the shock, and slowly recombines from the hydrogenic ion to the helium-like one. In this case the 
forbidden line of the triplet will be pumped by the recombination. A resolving power E/ AE > 60 
would be needed to separate these lines. The hydrogenic ion, however, also has emission lines- 
specifically, the Lya lines, which are at 1.47 keV, 2.01 keV, and 2.62 keV respectively for Mg Xii, Si 
XIV, and S xvi. In the recombining plasma, these lines may be an order of magnitude stronger than 
at equilibrium. Energy resolution of only E/ AE > 11 is required to distinguish them. The Einstein 
SSS has a FWHM of ~ 0.14 keV (Giaconi et al. 1979) which would allow at best a marginal hmit 
be placed on these lines. 



Table 5: The observed and calculated emission in the soft and hard ROSAT bands. 



Band 


Observations 


Model A 


Model B 


Model C 


soft (R4+R5+R6) 


2.026 ±0.069 


2.22 


0.923 


2.37 


hard (R7) 


1.838 ±0.067 


2.12 


0.915 


2.28 



Total Count Rate 

Table I shows the total count rates in the two ROSAT PSPC bands, soft (R46) and hard (R7), 
for the observations and the three different models. The observational values are extracted via the 
method of Snowden et al. (1996), while the model values are calculated using a response matrix 
generated through FTOOLS. The assumed Nh was 1.89 x 10^^ cm~^ for models A, B, and C. Using 
Xselect (Turner, 1996) and XSPEC (Arnaud, 1996) to reduce the data gives an ~ 10% higher count 
rate for these bands than the Snowden et al. (1996) method. The reason for this discrepancy is 
not completely understood, but is probably due to the improved background subtraction done in 
the Snowden et al. (1996) method. 
Reduction of Thermal Conductivity 

As discussed in Section 4.4, there is some controversy regarding the efficacy of thermal con- 
duction in the presence of a magnetic field. In order to explore the efi'ect of lowering the thermal 
conductivity, we calculated Model E, identical to Model A except with the thermal conduction 
constant [3 reduced to 10% of the standard value of 6 x 10~^ cgs. This model has a higher cen- 
tral temperature (Tc — 9.8 x 10® K), lower central electron density (uc — 0.51 cm~^), and steeper 
gradients than does case A. As a result, the emission spectrum (Figure is harder. The final fit 
is poor, with a reduced = 3.4. The PSPC model is overluminous by 5% while the SSS model 
is undcrluminous by a factor of 2.13. The best-fit value has — 1-78 x lO^^cm"^, slightly less 
than for the best fit to Model A. 



7 Discussion 



The methodology of this investigation has been the application of analytic theory as far as possible, 
followed up with numerical hydrodynamic modeling to verify the details. The main successes are 
summarized in the abstract, while the overview of the paper's organization and an index to the 
various results are provided in the introduction. A number of analytical tools, some in review 
and some unique and new, were presented in The synchrotron flux available from the van der 
Laan mechanism is discussed in §3.4.1, with the main result in Equation |2^. Proof that thermal 
conduction should lead to a fixed ratio of central to a mbie nt density at the cooling epoch, and the 
dependence of that ratio on conductivity appear in §3.3. The theoretical section is summarized 
in § |3.9| , and a schematic of the basic model appears in Figure The remainder of this section 
provides a brief summary as a backdrop for a fuller confrontation with a few key points. Each 
subsection presents successes first, followed by caveats, and ideas for observations. 



7.1 Energy, Environment, Age 

The model we have presented for W44 is the remnant of a 10^"'^ erg supernova. It occurred in a 
region with an explosion site density of about 6 cm~'^ and a density gradient whose scale height 
equals one remnant diameter at the present time. Our hydrodynamic studies assumed a smooth 
density distribution, but there must be irregularities to achieve the shell corrugation implied by 
the filamentary radio continuum emission pattern. The age of the remnant is that of the associated 
pulsar, about 20,000 years. 

None of observations we interpreted required the remnant to be interacting with dense molecular 
material in the vicinity. The model's dense end has a preshock density of only 10 cm~'^. The 
remnant itself shows no significant morphological relationship between its radio continuum and the 
oft mentioned Wooten cloud. It seems to be in a region with surrounding molecular clumps, but is 
not definitely encountering any. That said, we stress two caveats. First, the shell corrugation is a 



consequence of a lack of smoothness of the ambient medium on 1 to 10 pc scales and we have not 
explored the degree of density contrast required to generate that corrugation. Second, we found 
that the gamma ray flux could be raised to the observed level from the vicinity if the cosmic rays 
in the dense shell were interacting with 2 to 4 times as much mass as we predict. It is not known 
that the gamma rays do derive from the remnant, but, if it can be shown that the eastern shell is 
the source and has the presently quoted flux, it may be necessary to revisit the possibility of dense 
inclusions or high contrast irregularities. 



7.2 The Hot Interior 



We assumed that thermal conduction in the hot interior occurs at its unqucnched rate, limited by 
saturation (§4.4), with the result that the minimum interior density drops only to about 13% of 
ambient and the central temperature is limited to about 6 x 10*^ K. These conditions produce the 
observed central x-ray emission (5:3.3, 5.2.1, 3.2.1). The remnant's temperature is low in the outer 
parts, in part because it has become radiative. As a result, its center filled x-ray morphology is not 



accompanied by a bright shell, even prior to absorption (§5.2.4, 6.2.2). After absorption it yields 
approximately the observed intensity and spectrum. 

At normal abundances, the x-ray emission from the remnant interior is calculated to arise 
from approximately 30 Mq . One expects that the enriched stellar ejecta should have distorted 
the spectrum, the surface brightness distribution, and the overall intensity. And yet, our basic 
model reproduces most of the spectral detail, missing only noticeably in the stren gth of the Si XIII 
complex around 1.85 keV, and provides essentially the correct total count rate (§ |6.2.3 ). 



The only hints that the ejecta play a role are in the surface brightness distribution (§ |6.2.2| ). 
That distribution has three characteristics of interest. Unlike that provided by our basic model, the 
radial distribution has a relatively sharp central peak and is concave upward elsewhere. In addition, 
the image is somewhat lumpy, a characteristic which may have contributed to the idea that there 
are discrete evaporating clouds associated with the emission. Some of this distribution may arise 
from variability in the absorption screen, but we are inclined for reasons given in the text to believe 
that at least the central peaking is real (§2.5). Our ad hoc experiment with moderately enhanced 
abundances in the stellar ejecta (§ |6.2.2 ) easily provided sufficient additional central emission, but 
far more can be done to clarify the situation. The observed distribution and spectrum depend 
at least on the abundance distribution in the stellar envelope, on the wind activity of the star in 
shaping its environment prior to the explosion, on the degree to which dust formation extracts 
elements from the gaseous phase, on the instabilities in the wind bubble and explosion activity 
tending naturally to form lumps, on the potential for some of those lumps to cool and return to 
high density, on the radial mixing profile of the interpenetrating fiuids of ejecta, wind and ISM, 
on the thermal conductivity within the stellar wind bubble and in the hot interior of the remnant, 
on the electron-ion thermal equilibration rate, on the dust destruction rates, and on the stellar 
density profile and the irregularities in the medium which will influence the focussing and location 
of the reverse or reflected reverse shock at the present time (§ |3.2.1 ). 

Improved imaging and spectra of W44 X-rays can very much clarify the const rainin g infor- 
mation, particularly on the ionization history and current temperature of the gas (§ |6.2.3 ), but it 
would be premature to suppose that existing modeling efforts have been sufficient for meaningful 
comparison with data taken from the central third of the remnant. A great deal more careful 
work lies ahead, with major difficulties and dangers in the sophistication needed in approximating 
the interrelated phenomena. The good news is that many models of the details will be totally 
unacceptable, producing far more X-rays with very much the wrong spectra. The payoff will be 
multiparameter constraints on the possible characteristics of the star, it's interaction with the sur- 
roundings, and on the degree of thermal conduction allowed or required. As for the concavity of the 
angle averaged surface brightness distribution, modelers must take care to be sure that scattering 
of the X-rays by intervening gas is included in their final comparisons. 



7.3 The Shock and Dense Shell 

The model remnant's shock shows cooling effects everywhere, but a dense shell has formed only on 
the higher density end (Figure d and §|t], I^, |5.2.1D . With that end tilted somewhat away from 



the line of sight, the half shell can be seen as a receding distribution of HI whose velocity structure 
agrees with observations (§5.2.2). The radio continuum arises from the dense shell, via synchrotron 
emission in an approximately 200 /iG field (§3.4.1, 5.2.3| ), with its filamentary structure due to local 
tangencies of the corrugated shell with the line of sight (§5.2.3). The shell should also be a source 
of gamma rays, from both bremsstrahlung by the CR electrons and 7r° production by the protons, 
at about l/4th to 1/2 the rate measured by EGRET as arising from the vicinity (§ |3.6| ). 

The radiative shock and ongoing shell formation generate the usual radiative remnant's spec- 
trum, in UV, opt ical, and IR, with rates that correspond to recent observations of Ha, [SII], and 
63 /im [Ol] (§3^). Much of this emission is associated with the edge of the dense shell and should 
correlate well with the filamentary structure, though some could be more diffusely arranged in the 
remnant's equatorial belt where shell formation is currently strongest. Both aspects are consistent 
with existing observations, though in the optical the observations are made through a variable 
screen of several magnitudes of foreground extinction (measured by the absorption column of the 
x-ray spectrum). 

Because the filaments arise geometrically from tangencies to the shell, any correlated line emis- 
sion should appear at the local rest velocity (expansion is perpendicular to the line of sight), a 
characteristic that distinguishes filaments due to tangency from those due to actual filamentary 
(rope- like) structures. This velocity null for the filaments has been confirmed in the associated 
OH maser emission (Frail, 1998, private communication), but should be checked in Ha because of 
possible disagreements over the origin of the maser emission. 

Unfortunately, the large and variable optical depth to W44 degrades the information available 
from optical intensity studies. But a velocity resolved study of the Ha, or of other spectral com- 
ponents potentially associated with the shell and or evaporating clouds, could be useful in several 
ways. We expect emission from the existing shell, from the immediate surroundings, and from the 
equatorial belt where shell formation is in progress. The shell component should strongly resemble 
the HI distribution in velocity and space, but can likely be detected from the approaching (shell 
forming) side as well as from the receding side, and traced to the shell edge at low velocity, the 
latter not being possible in 21 cm due to confusion with other material along the line of sight. In 
addition to the important verification of the the zero radial velocity of the filamentary emission, a 
search is needed for a possible emission component associated with the x-ray bright spots, as would 
be expected from evaporating (or cooling and condensing) lumps. Perhaps most interesting is the 
possibility of seeing emission from gas which is just joining the inner surface of the shell. If the 
cooling is sensitive and therefore somewhat erratic, the condensation flow could be irregular and 
have some components expanding at a noticeably higher speed than the shell they are overtaking. 

We found that the van der Laan mechanism, sweeping up of ambient cosmic rays and magnetic 
field and compression of both into the dense shell, can p lausibly provide both the observed spectrum 
and intensity of the radio synchrotron emission (§ [3.4. l| ), but there was a hidden difficulty with the 
formulation that amounted to assuming that the ambient CR electron density there is about 4 times 
greater than that estimated for the solar vicinity. This is not necessarily a problem. Attempts 
to measure the synchrotron emissivity of the general ISM have been difficult and seem to have 
both significant uncertainty and genuine variability. In addition, W44 is in an active region at 
the base of the Aquila supershell (Maciejewski et ai, 1996) and could very well be encountering 
an enhanced CR population. Other important processes potentially contributing to the cosmic 



ray electron population are discussed in §3.4.2. A caveat in our analysis is that our use of the 
viewing-angle averaged synchrotron radiation formula may have been misleading in the presence 
of the orderly tangential magnetic field seen in W44. 

The extremely high angular resolution of the OH maser locations and their relationship to radio 
continuum ridges and bright spots (Claussen et al. 1997) appears to us to require the masers to arise 
from the leading edge of the dense shell, i.e. at the trailing edge of the cooling and recombination 
zone of the radiative shock, with masering occuring tangential to the shell's outer surface so that 
it is observed only along the filaments. Consistently, the magnetic field measured via the maser 
emission is just that found for the shell in our model (§3.4.2), and as noted above, their radial 
velocities are at the local rest velocity. We are not aware of models for shocks in such low density 
gas to test for the presence of sufficient OH, or of maser pumping models for this nonequilibrium 
environment, but we note that the models by HoUenbach and McKee (1989) for shocks in denser 
gas do have a significant abundance of OH on leaving the photoionization/rccombination plateau. 



and show a local peak in the OH emission at a column density behind the shock of about 10^*^ 
cm^^, where the temperature is dropping rapidly through 1000 K. This result may not be sensitive 
to the initial density (McKee, private communication 1998); and the required column density is 
achieved in our model. The anticipated neutral density in this region is a few hundred cm~'^, 
limited by the magnetic pressure. The observations cry for studies of 100-150 km s~^ shocks into 
densities of 10-100 cm~'^ and their potential for creating molecular gas and generating OH maser 
activity tangential to the shock surface, possibly with the pumping arising in transience of the 
material (e.g. formation of the OH in excited states). The solution could be elusive, however, 
based on the observation that there are many small maser spots but far more filamentary regions 
without masering, at least not beaming towards us. 

Our predicted shell mass would also agree better with the HI data if much of the shell were 



molecular (§ 5.2.2 ). If there is a strong molecular line which is insensitive to reddening and asso- 
ciated with radiative shocks, perhaps from the fluorescence region at the outer edge of the dense 
shell, a line not present in the ambient medium or foreground molecular gas, it could provide an 
excellent probe of the shell's actual extent and velocity structure. The model's overall shell mass 
is about 450 Mq of which perhaps half is molecular. If it can be shown that the molecular mass in 
the shell is substantially larger than this, perhaps in structures related to the shell corrugation, it 
could also remove the apparent discrepancy between our model's gamma ray production and the 
EGRET measurements, while requiring only some adjustment of the ambient density distribution 
in the next generation of models. 

7.4 Distance and Fine Tuning 



The oft-quoted distance for W44 is 3 kpc, but we reviewed this assessment in §2.2, concluding that 
the appropriate modern number is 2.5 to 2.6 kpc. Our theoretical attempts to choose parameters 
that would produce the correct angular size at this distance failed slightly, forcing us to describe the 
model characteristics for a distance of 2.2 kpc. It is not impossible that the remnant is this close, and 
within the uncertainties of such modeling, the differences are inconsequential. Nevertheless, it may 



be worth comment. Section 3.S introduced scaling laws for adjusting our results to remnants located 
at different distances. Two scalings are presented, exact Sgro-Chevalier scaling, and approximate 
scaling for constant explosion energy. For Sgro-Chevalier scaling to 2.5 kpc, we found that the 
only appreciable changes were a 12% reduction in the x-ray surface brightness and luminosity, and 
a factor of 1.3 increase in explosion energy and shell mass. The increase in explosion energy is 
the only uncomfortable part, leading us to favor the more approximate scaling. In that case, the 
only appreciable change is the reduction of the x-ray surface brightness and luminosity by a factor 
0.53. Recalling that Model A was about 25% overluminous leads to no worse overall agreement. 
However, considering that a significant part of the total x-ray emission may be associated with 
the elementally enriched center as discussed above, this increase in distance is an improvement, 
allowing a greater degree of enrichment of the interior without having to resort to a reduced thermal 
conductivity for compensation. With this model, the explosion site density would be only 4.6 cm~'^ 
and the age 23,700 years. Both arc intermediate between Model A and the parameters found by 
Harms, et al. (1997). 

7.5 Relationship of the Model to Other Remnants 

Finally we would like to remark on the possibility of generalizing this model to other remnants. 
To have both a strong radio shell component and an interior which is luminous in thermal X- 
rays requires both a substantial compression of the dense shell, and a substantial interior density 



and pressure. In §3.4.1 we found that the radio synchrotron surface brightness shortly after shell 
formation should be roughly proportional to Kpsheii , while in § p.8| we found the x-ray surface 
brightness of the interior proportional to Rpf^jg^j^,^. We have not explored our hydrodynamic 
models to learn how rapidly these surface brightnesses fade as the remnant continues to expand, 
but we expect a very considerable drop in the x-ray emission by the time it has grown very far. The 
radio continuum is less sensitive in the first place, and will drop more gradually because it rises with 
cooled mass and because remnants in inhomogeneous media have shell formation occuring over part 
of their surface for a protracted period of time. We naturally expect then that remnants with radio 



shells will be more common than those with x-ray luminous interiors, and can concentrate on the 
conditions needed for the latter. At shell formation, when the outer edge is too cool to radiate in 
x-rays, the surface brightness of the interior should be roughly proportional to RJoq/- In that case 
it is clear that the class of remnants to which our model applies should have a very steep dropoff 
of surface brightness with size at shell formation, as well as with size after shell formation. Using 
§3.2 to relate shell formation radius to ambient density, we infer that the x-ray surface brightness 
of the interior at shell formation is approximately proportional to n^. As we have not explored the 
effects of any presupernova wind on the later internal structure of the remnant, this result may 
be somewhat naive, but our tentative expectation is that an x-ray bright interior for a supernova 
remnant of typical energy will occur only when the explosion occurs in an extended region of higher 
than average interstellar density. 

For contrast, consider the Cygnus Loop. Its radio continuum derives largely from compressed 
material behind roughly 100 km s^^ radiative shocks moving into recently encountered large scale 
clouds of density roughly 8 cm^"^, not so different from W44. Its flux is similar to W44's, its 
distance about a third as great, so its radio continuum luminosity is lower by about a factor of 10. 
The remnant volume is about 4 times that of W44, so the mean internal pressure is about 4 times 
lower. The brightest x-rays are associated with the interaction with the large dense clouds, and 
are found just inside the radiative shocks. The distribution of the more diffuse x-rays, however, 
resembles that expected from a Sedov structure with ambient density 0.16 cm~^. We therefore 
estimate the interior density to be about l/30th that of W44 and the interior temperature about 
7.5 times greater, 4 to 5 xlO^ K. Continuing to assume that the abundances are unaffected by the 
ejecta, which provides the correct count rate for W44, and consulting our ROSAT PSPC counts 
per emission measure tables with negligible absorption for the Cygnus Loop and the amount found 
above for W44, we conclude that the observable (after absorption) surface brightness of the Loop 
interior should be lower by roughly a factor of 10. That reduction appears modest until compared 
with the radiation from the 30 times denser edge. At a temperature of 2.4 xlO^ K, the "shell" 
emission is about 2 orders of magnitude brighter in projection across the center of the remnant 
and much brighter yet as seen near the edge. Although there are radiative shocks on part of the 
surface of the Cygnus Loop, much of the evolution is in low density and is very far from having the 
low edge temperatures associated with incipient cooling and shell formation. Only an instrument 
with sensitivity to very hot gas could discern the high temperature interior from the softer "shell" 
emission. 

There are two ways of creating much larger remnants with x-ray emitting interiors. One, ex- 
plored by Shelton (1998) in the context of halo and interarm supernovae, simply involves exploding 
in a very low density (0.01 cm^^) and waiting for the edge temperature to drop too low to show up 
as an x-ray shell. The other was explored in Smith and Cox (1998) in connection with modeling the 
Local Bubble. A second supernova within the cavity of an earlier remnant can reheat the interior 
and bring it back into x-ray emission. And that emission can last a substantial length of time, 
possibly making unlikely events yield not uncommon objects in the sky. The two approaches are 
not entirely dissimilar in that the Local Bubble model starts with a higher ambient density and 
uses the first explosion to create the low density into which the second explosion initially evolves. 
One difference is that the multiple explosion model has denser surroundings which slow the shock 
and lead more quickly to a low temperature edge. 



7.6 Closing Remarks 

This study was undertaken to show that the characteristics of a heavily observed and very excit- 
ing supernova remnant, W44, could be understood in the context of a very simple model, a model 
which did not invoke evaporating clouds in the center, or blast wave accelerated clouds to make the 
observed HI, or interaction with a dense molecular cloud environment, or diffusive shock injection 
and acceleration of cosmic rays. The required features were a relatively dense ambient environ- 
ment, a significant density gradient, shell corrugation due to moderate irregularities in the ambient 
medium, betatron acceleration of the ambient cosmic ray population, and a significant amount of 
thermal conduction in the hot interior. It was hoped that the observed intensity of x-rays from 
the hot interior would then finally give us a relatively conclusive measure of the magnitude of the 
thermal conductivity. 



The model was wildly successful, but not without its disappointments. We were able to show 
that compression of the ambient cosmic rays could lead to approximately the correct radio contin- 
uum flux and spectrum, and to a signiflcant contribution to the gamma ray intensity, but could 
not show that there is no room left for diffusive or other forms of acceleration. We found that 
the intensity and spectrum of the interior x-ray emission were consistent with unquenched thermal 
conduction, but could not convince ourselves that the metallicity contamination by the cjccta and 
preexplosion evolution could be neglected. As a result we are less confident in our determination 
of the magnitude of the thermal conductivity in the hot interior of the remnant. On the positive 
side, the need for clouds and clumps in the remnant interior vanished, as did any hint of interaction 
with dense molecular clouds. The shell emission characteristics were in beautiful agreement with 
the measurements, and comparison with the OH maser observations led to the ideas that the shell 
could be largely molecular and that the trailing edge of the cooling region is the site of the masering 
activity. 
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10 Figure Captions 



Figure 1: Schematic of the Model for W44. Shows annotated section through remnant along plane 
containing the look direction. Inset shows projection of half shell on the sky. 

Figure 2: Number density of each cell in the 2d model, projected on the symmetry axis. The 
ambient density distribution, partial shell formation, factor of 4 compression at the tenuous end's 
shock, and the moderately uniform and non- negligible central density are all apparent, as are the 
regions of finer gridding in the adaptive mesh. The remnant age is 20,000 years. 

Figure 3: Surface plot of the 2d model's temperature distribution at 20,000 years. The temperature 
structure is approximately parabolic, fairly uniform in the central half (in radius) with a steady 
drop in the outer half to near zero at the edge. On the tenuous end, the post shock temperature 
of the shock can be seen, along with the trough behind it showing its transition to radiative. 

Figure 4: Surface plot of the 2d model's pressure distribution at 20,000 years. The constant 
central pressure, sloping shoulders in the outer half, strong variations in edge pressure over the 
surface and the deep trough of cooled but so far uncompressed material can be seen. (The spikes 
are an artifact of the figure preparation, not a characteristic of the model.) 

Figure 5: (top-right) HI mass from the simulated SNR viewed at 40° from the symmetry axis. 
The axis is also tilted toward the southwest, (bottom-left) The shaded contours represent the HI 
brightness observed for W44 between wlsr = 135 and 170 km s"^ (Koo and Heiles 1995). The dots 
represent HI mass from the model for the same velocity interval, assuming a 43 km s~^ systemic 
velocity, (^top-left) Position velocity diagram for the vertical aperture indicated in the top right 
panel, (bottom-right) Position velocity diagram for the horizontal aperture indicated in the top 
right panel. 

Figure 6: (top) Position-velocity diagram of the fast, receding HI gas. The shaded contours 
represent HI brightness for a cut along aigso — 18'^53™30^, from Koo & Heiles (1995). The dots 
represent HI mass for the vertical cut indicated in the previous figure. A systemic velocity of 43 
km is assumed, (bottom) Same as top, for a horizontal cut through Sig^o = -|-01°18'. 

Figure 7: Map of the synchrotron emission from the simulated SNR. In the upper right panel, the 
remnant is viewed at the same orientation as in Figure ^ whereas the lower left panel shows a view 
perpendicular to the symmetry axis. The upper left and lower right panels show the integrated 
intensity profiles for the horizontal and vertical slices shown. 

Figure 8: Map of the synchrotron emission from the simulated SNR, in which a pattern of radial 
crinkling has been superposed to illustrate the formation of apparent filaments. The remnant is 
viewed at the same orientations as in Figure |^. 

Figure 9: The 2d model x-ray surface brightness distribution. Includes both hard and soft ROSAT 
PSPC bands, at a column density of 2x10^^ cm~^. 

Figure 10: The 2d model x-ray surface brightness distribution with no absorption. The 2d model's 
thermal X-rays are not edge brightened (shell like) even without absorption. 

Figure 11: Density distributions of Id models at 20,000 years. With thermal conduction (Model 
A, solid line), the electron density is almost 1 cm^'^ in the center and changes slowly with radius 
except near the edge. Without thermal conduction (Model B, dashed line) the density is extremely 
small at the center and increases rapidly with radius. The structures near 11 pc show the early 
stages of shell formation. 



Figure 12: Temperature distributions of Id models at 20,000 years. With thermal conduction 
(Model A, solid line), the temperature is 6 x 10^ K in the center and the changes only slowly with 
radius except near the edge. Without thermal conduction (Model B, dashed line) the temperature 
is very high in the center and changes rapidly with radius. The bump at 6 pc is an artifact of 
the code associated with the contact discontinuity. The structures near 11 pc show the early stages 
of shell formation. 



Figure 13: Thermal pressure distribution of Id models at 20,000 years, with (Model A, solid line) 
and without (Model B, dashed line) thermal conduction. 



Figure 14: a: The ROSAT PSPC Soft Band (R4+R5+R6, 0.4 - 1.6 keV) count rate for Model A 
(with thermal conduction, solid line). Model B (without thermal conduction, dashed line), and the 
observed rates toward the north (cross), east (star), south (diamond), and west (triangle). Models 
A and B were attenuated with an absoption column density of Nh = 1.89 x lO^^cm"^. b: Same as 
a, but for the ROSAT PSPC Hard Band (R7, 1.0 - 2.0 keV). 



Figure 15: ROSAT PSPC (a) and Einstein SSS (b) spectra fit with Model A. 



Figure 16: ROSAT PSPC (a) and Einstein SSS (b) spectra fit with Model C. 



Figure 17: ROSAT PSPC (a) and Einstein SSS (b) views of Model A comparing coUisional ioniza- 
tion equilibrium (dashed) and full non-equilibrium ionization (solid). The non-equilibrium model 
is the same as shown in Figure |l^. 



Figure 18: ROSAT PSPC (a) and Einstein SSS (b) spectra fit with Model E, identical to Model A 
but with the thermal conductivity reduced by 10. 
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